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CHAPTER 1 


Background and results 


In recent years, substantial progress on a statement by Vladimir Igorevich 
Arnold concerning the stability of the planetary system has been achieved [22, 
2, 23, 33, 20, 14, 27, 9], 

It sounds as follows. 

“For the majority of initial conditions under which the instantaneous orbits of the 
planets are close to circles lying in a single plane, perturbation of the planets on 
one another produces, in the course of an infinite interval of time, little change 
on these orbits provided the masses of the planets are sufficiently small. [... ] In 
particular [... ] in the n-body problem there exists a set of initial conditions having 
a positive Lebesgue measure and such that, if the initial positions and velocities 
of the bodies belong to this set, the distances of the bodies from each other will 
remain perpetually bounded.” [2, Chapter III, p. 125]. 

Solving the differential equations of the motions of the planetary problem, i.e., n 
planets interacting among themselves and with a star via gravity is, for n > 2, 
a problem with ancient roots. This story goes back to Sir Isaac Newton - who 
brilliantly solved the case of two bodies and then, tackling the analogous one for 
three bodies, soon realized the necessity of turning to a “perturbative” study (except 
for naming it a “head ache problem”) - passed through investigations by eminent 
mathematicians like Delaunay, Lagrange, the prize publicly announced by King 
Oscar II of Sweden and Norway and awarded to Henri Poincare, but its “solution” 
is nowadays open. Chaotic and stable regions may coexist [2, 17, 11]. 

The question received a new mathematical description, and a strong modern en¬ 
dorsement, after A. N. Kolmogorov announced, at the International Congress of 
Mathematicians of 1954 in Amsterdam, what is now almost unanimously considered 
the most important result of the last century for dynamical systems: the theorem 
of conservation of the invariant torus. This breakthrough result, next enriched by 
substantial contributions by J. Moser and V. I. Arnold himself [22, 26, 1], states 
that for a generic Hamiltonian system close to an integrable one, i.e., a system of 
the form 

=h{I) + fif{I,g>) {I,g>)eBxT^ BCR^ T:=R/(27rZ) p < 1, 

the major part of unperturbed motions survives, after a small perturbation is 
switched on, provided suitable “non-degeneracy” conditions are verified by the “un¬ 
perturbed part” h. Moreover, the theory provides precise arithmetic (“diophantine”) 
properties to be verihed by the “unperturbed frequencies” w* = i9h(/*), so that they 
will be preserved in the full system. 

In 1962, V. I. Arnold, extending Kolmogorov’s ideas, and looking for an applica¬ 
tion to the planetary problem, at the International Congress of Mathematicians of 
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Stockholm, announced the theorem of stability of planetary motions quoted above. 
In 1965 Kolmogorov and Arnold were awarded the Lenin Prize for their studies on 
the stability of the planetary problem - but the story was not finished there. 

In order to introduce the results of this paper, we highlight basic facts of this story 
and its continuation, referring the reader to [16, 5, 28, 10, 29] for more notices. 

The planetary problem is close to the integrable problem of n uncoupled two-body 
problems, where each planet interacts separately with the sun. The mutual interac¬ 
tions among planets are regarded as a perturbing function, the smallness of which 
is ruled by the planets’ masses. However, as a perturbed system, the planetary 
problem has a limiting degeneracy. Its associated integrable system (the two-body 
problem) is “super-integrable”: it has more integrals than degrees of freedom. At 
a technical level, the limiting degeneracy is exhibited by the disappearance of de¬ 
grees of freedom in the unperturbed part. Therefore, continuing the unperturbed 
motions to a positive measure set of quasi-periodic trajectories might, in general, 
be not possible, in absence of further informations on the perturbing function. 

Arnold found, for the planetary problem, a brilliant solution to the problem of the 
limiting degeneracy. This led him to add to the assumptions and assertions that 
are proper of perturbation theories (e.g., “the masses of the planets are sufficiently 
small”, “set of initial conditions having a positive Lebesgue measure”, “the distances 
.. .will remain perpetually bounded”)a further requirement of smallness of eccen¬ 
tricities and inclinations of the unperturbed Keplerian ellipses (“the instantaneous 
orbits of the planets are close to circles lying in a single plane”). Let’s summarize 
Arnold’s ideas. 

Choosing, as Arnold did, Poincare coordinates [30] (see, also [2, Ch. Ill, §2], or, 
e.g., [8, 15]), the system takes the usual close to be integrable form 

Up02 — fiKep T fi fpoi: 

where /i is a small parameter related to the planetary masses, but the unperturbed 
“Keplerian” part hKep(A) depends on only n action variables A = (Ai, • • • , A„) (re¬ 
lated to the semi-major axes of the instantaneous Keplerian ellipses), out of an over¬ 
all of 3n degrees of freedom. The perturbing function, fpoi, on the other hand, de¬ 
pends on all the coordinates: the actions A, their conjugated angles £ = (£ 1 , • • • , in) 
(proportional to the areas of the elliptic sectors spanned by the planets), and, more¬ 
over, on some other coordinates (p, q) = (pi, • • • , p 2 rn qi, • ‘ ‘ j ^2n), 4n-dimensional, 
related to those (“secular”) quantities (eccentricities, inclinations, nodes and peri¬ 
helia of the ellipses) that in the unperturbed problem stay fixed, and for this reason 
do not appear in hxep. 

It is of great help that the averaged perturbing function (with respect to the angles 
^) /po 2 (A, p, q) enjoys several parities in the coordinates (p, q), geometrically related 
to its invariance by rotations and reflections with respect to the coordinate planes. 
The “secular origin” (p, q) = 0, corresponding to all the planets moving on ec¬ 
centric circles in the same plane, turns out to be an elliptic equilibrium point for 
the averaged perturbing function, for any value of A. 

Arnold brilliantly argued to exploit this circumstance to his purpose. By Birkhoff 
theory, one might think to switch to another set of canonical coordinates (A, £, p, q), 
analogous to Poincare’s coordinates, possibly defined only for (p, q) in a small neigh¬ 
borhood of radius e around the origin, such that the Hamiltonian of the system, or. 
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more precisely, its ^-averaged (“secular”) perturbing function /Bir, takes a “normal¬ 
ized form” : it is a polynomial, fBir,tr, of some degree greater or equal than two in 

~2 I —2 

the combinations (“degenerate actions”) ^ i = 1, • • •, 2n, plus a remainder 

with a higher order. Roughly, Arnold projected to solve the limiting degeneracy 
by conjugating the planetary system to a new system, whose unperturbed part was 
just the truncated, normalized Hamiltonian 

hKep “k fnir^tr 

SO as to recover the standard set up of KAM theory. With these ideas in mind, he 
proved the following impressive result and next applied it to the planar three-body 
problem. It states that stable trajectories occupy a positive measure set of the phase 
space, and are more and more dense closer to the elliptic equilibrium. Hence, the 
smaller eccentricities and inclinations are, the larger the number of stable motions 
is. 

‘The Fundamental Theorem” (V. I. Arnold, [2]) If the Hessian matrix of h 
and the matrix of the coefficients of the second-order term in Ti in fsir (“torsion”, 
or “second-order Birkhoff invariants”) do not vanish identically, and if fj, is suitably 
small with respect to s, the system affords a positive measure set of quasi- 

periodic motions in phase space such that its density goes to 1 as e —>■ 0. 

Arnold perfectly knew that, in order to apply the Fundamental Theorem to the 
problem in space, one should previously treat an unpleasant fact: one of the first 
order Birkhoff invariants vanishes identically. He was aware that the reason for this 
first-order degeneracy was to be sought into the existence of two non-commuting 
integrals, the two horizontal components of the total angular momentum of the 
systems. If, apparently, a vanishing eigenvalue strongly violates the construction 
of the normalized system (a deeper analysis of the symmetries of the perturbing 
function [25, 8], however, shows that the identically vanishing eigenvalue is not a 
real obstruction), a major problem definitely prevents the application of the Fun¬ 
damental Theorem: an infinite number of coefficients of any order of the (formal) 
Birkhoff series vanishes identically, among which one entire row and a column in 
the torsion matrix, which so is identically singular, and the reason is again the in¬ 
variance by rotations. The proof of this generalized degeneracy is in [8]. We recall 
here that even Herman had raised a question about the degeneracy of torsion [20, 
p. 24]. 

We do not know weather Arnold was aware of the infinite degeneracy of the nor¬ 
malized system (he did not even mention the vanishing of torsion in his paper). 
He however suggested two different strategies for the three- and the maiiy-body 
case, of which he provided very few and somewhat controversial details. As for the 
three-body problem (his ideas for the many-body case will be recalled a few below), 
he proposed to reduce the integrals (hence, the number of degrees of freedom) of 
the system by switching to a system of canonical coordinates going back to the XIX 
century, worked out by Jacobi and Radau [21, 32], which in literature go under the 
name of Jacobi reduction of the nodes. The idea was later completely developed by 
P. Robutel [33], who, in a deeply quantitative study, checked the non-degeneracy 
assumptions required by the Fundamental Theorem. 

Finding a system of canonical coordinates that do the job of Jacobi reduction of the 
nodes when the number of bodies is more than three has been a central difficulty 
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for a long time [2, 25]. At this respect, Arnold sadly commented: “In the case of 
more than three bodies there is no such elegant method [as Jacobi reduction of the 
nodes] of reducing the number of degrees of freedom.” ]2, Ch. Ill, §5.5, p. 141]. 

Exactly twenty years later, F. Boigey and A. Deprit refuted this sentence ]3, 12]. 
They indeed were able to extend Jacobi-Radau reduction to the four-body and 
general problem, respectively. It should be remarked, anyway, that, while the 
works by Jacobi, Radau and Boigey provide canonical coordinates on suitable sub¬ 
manifolds of the phase space, the one by Deprit is more general and clarifying, since 
it provides a set of canonical coordinates for the whole phase space and allows us 
to recover his predecessors by restriction. 

The utility of Boigey-Deprit’s coordinates was not suddenly clear. Neither Boigey 
nor Deprit ever provided any motivation of their study, or foresaw applications. 
The only application that is known to the author up to 2008, concerning indeed 
Deprit’s coordinates, stands in a paper by Ferrer and Osacar, in the 90s, to the three 
body problem ]18]. But this case is not really exhaustive, since for three bodies 
Deprit’s and Jacobi-Radau’s coordinates coincide. A reason why Boigey-Deprit’s 
coordinates have been forgotten so long might be that, for more than three bodies, 
they actually have a less natural aspect, compared to the classical case of Jacobi. A 
sort of “hierarchical” structure in the geometry of Deprit’s coordinates discouraged 
the author himself, who, at the end of his paper, declared: “Whether the new phase 
variables are practical in the general theory of perturbation is an open question. At 
least, for planetary theories, the answer is likely to be in the negative. But finding 
a natural system of coordinates for eliminating the nodes in a planetary cluster was 
not the intention of this note.” [12, p. 194]. 

In the meantime, in 2004, the first general proof of Arnold’s stability statement ap¬ 
peared. It was by Jacques Fejoz, who completed investigations by the late Michael 
Herman ]14] - but the different procedure that Herman had in mind did not rely on 
the necessity of handling, explicitly, good coordinates. Indeed, Herman conceived a 
proof based, besides on a “twist-less” KAM theory going back to H. Russmann ]34], 
on indirect arguments of Lagrangian intersections in order to bypass the so-called 
“secular resonances”. See [10] for more details. 

In 2008, Boigey-Deprit’s coordinates were rediscovered by the author ]27], in a 
slightly different, “planetary” form. The rediscovery was motivated by the purpose 
of realizing Arnold’s program (i.e., applying the Fundamental Theorem quoted 
above directly to the planetary Hamiltonian) in the general case, so as to obtain a 
detailed information about the tori frequencies, the measure of the invariant set and 
the symplectic structure of the phase space. The utility of Boigey-Deprit’s coordi¬ 
nates became suddenly clear: switching (in order to overcome certain singularities of 
the chart) to a regularized version, called “RPS” coordinates, (acronym standing for 
“Regular, Planetary and Symplectic”), allowed them to derive the Birkhoff normal 
form of the planetary problem, to prove its non-degeneracy, and hence to complete 
the application of the Fundamental Theorem to the general problem. These results 
have been published in ]6, 7, 9]. 

Qualitatively, RPS coordinates are very different from JRBD (Jacobi-Radau-Boigey- 
Deprit); rather, they are more similar to Poincare coordinates. The mentioned 
parities and the elliptic equilibrium of the averaged system are still present in the 
RPS-averaged system. But, as an advantage with respect to Poincare coordinates. 
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the RPS perform^ a “partial reduction” of the rotation symmetry - in contrast 
with JRBD coordinates, which reduce “fully”. This way, all the degeneracies of 
the Birkhoff series mentioned above are removed at once, and the non-degeneracy 
assumptions of the Fundamental Theorem may be checked. 

We like to recall now Arnold’s strategy for the many-body case: more than forty 
years earlier, he foresaw to construct a system of coordinates analogous to RPS, 
via a Taylor series in Poincare coordinates [2, Ch III, §5, n. 5, p. 141]. 

Indeed, both the reduction of the nodes and this latter reduction are available 
whatever the number of bodies is. 

The possibility of switching from Delaunay-Poincare to the more fruitful JRBD, 
or even RPS coordinates, is an effect of the limiting degeneracy. This gives in 
fact the opportunity of remixing coordinates related to secular quantities, and, 
simultaneously, keeping the Keplerian term hxep unvaried. 

Following this idea, in this paper, we show that other systems of coordinates may be 
determined for the planetary problem which, as well as JRBD and RPS coordinates, 
are well adapted to overcome the degeneracy due to rotations, and, moreover, enjoy 
some different properties. 

We present a full reduction, which we call P-map, or perihelia reduction. It refines 
JRBD coordinates in two respects. 

Firstly, the P-map is well defined in the case of the planar problem, while JRBD 
coordinates are not. Everyone knows, in fact, that the starting point for the Radau- 
Jacobi reduction is the so-called “line of the nodes”, the straight line determined by 
the intersection between the planes of the two orbits. When the orbits of the two 
planets belong to the same plane, this is not defined. A similar circumstance arises 
for Boigey-Deprit’s coordinates, since their construction relies on certain straight 
lines in the space, which again lose their meaning in case of co-planarity. 

The proof of Arnold’s theorem given in [27, 9] is not affected by such singularity, 
since, as said, it relies on RPS coordinates, which, at the expense of one more degree 
of freedom, are well defined for co-planar motions - in that case they reduce to the 
classical Poincare coordinates. 

It has its consequences when one wants to compare results for the fully reduced 
systems, in space or in the plane. The singularity of the chart does not allow one 
to state that motions in the spatial problem with minimum number of independent 
frequencies starting with very small inclinations stay close to the corresponding 
planar motions. Notwithstanding further studies appearing in [28], where this 
problem is partially overcome (via the construction of regular coordinates for co- 
planar motions defined locally), it would be nice, in principle, to handle a global 
system of action-angle coordinates which completely reduces rotations and is shared 
simultaneously by the planar and the spatial problem. 

Secondly, the P-map is well adapted to reflection symmetries of the problem, while 
JRBD coordinates are not, as discussed in [25, 29]. 


^In the framework of the study of canonical coordinates for the planetary system, by “partial 
reduction”, we mean a system of canonical coordinates where a couple of conjugated coordinates 
consists of integrals (e.g., functions of the three components of the total angular momentum). 
By “full reduction”, we mean a partial reduction where also another integral appears among the 
coordinates. The terms “partial reduction”, “full reduction” have been coined in [25]. 
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Reflection symmetries are parities of the Hamiltonian expressed in Cartesian coor¬ 
dinates. As known, this does not change under arbitrary changes of the signs of 
positions or momenta coordinates. They are not related to integrals. Therefore, 
it might be a nice fact, and in general useful for applications, to have a system 
of coordinates that, after integrals are reduced, parities associated to reflections 
are maintained. Quite often parities are associated to equilibria, and equilibria to 
stable motions; an example is provided a few lines below. 

We shall apply the P-map by proving a variant of Arnold’s stability theorem. We 
shall face up to a question raised again by Arnold in his fantastic paper on the 
possibility of removing the constraint on eccentricities and inclinations. He indeed 
proved that, at least for the planar three-body problem, there is no need to assume 
their smallness. Rather, it is sufficient that the trajectories of the planets are away 
enough so as to avoid collisions. He obtained this stronger result by exploiting 
the convergence of the Birkhoff series associated to the averaged perturbation, a 
very particular and happy circumstance, due to the few degrees of freedom of the 
problem. 

From the mathematical point of view, the question is whetherstrategies exist for 
finding stable motions other than the one of exploring the neighborhood of the 
elliptic equilibrium. 

Concerning instead the physical relevance, asteroids or some trans-Neptunian ob¬ 
jects have motions with relatively large eccentricities and inclinations and an almost 
continuous spectrum of frequencies. 

Besides the mentioned stronger result by Arnold, some other statements in the 
same direction have been obtained for the case of the spatial three-body problem 
and the planar problem with any number of bodies [28]. Here, the measure of the 
invariant set has been estimated to be larger and larger as the planetary masses 
and the semi-axes ratios are small, no matter the smallness of the eccentricities 
and inclinations - the proof relies on an argument of convergence of a significant 
approximation of the Birkhoff series. Other results in this direction have been 
announced by J. Fejoz, since late 2013 [13]. 

Even though the arguments of ]2, 28] do not apply to the general spatial problem, 
since no significant approximation of the Birkhoff series associated to the averaged 
perturbation is integrable, using the P-map, we shall prove the following. 


Theorem A Fix numbers 0 < < 0.6627 .i — 1, ■■■ ,n. There exists a 

number N depending only on n and a number oq depending on e^, Ci, and n such 
that, if a < Oq, in a domain of planetary motions where the semi-major 

axes ai < 02 <■■■< a„ are spaced as follows 


{*) 


Rj < ai < af 


with 


a 


± 

1 


Q ,^( 2 "+ i - 2'*->+2 + l - i ) 


there exists a positive measure set the density of which in phase space can be 

bounded below as 

dens(K^,„) > 1 - (log 

consisting of quasi-periodic motions with 3n — 2 frequencies where the planets ’ ec¬ 
centricities Ci verify 

e.i < Ci < Ci. 
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Before we switch to details, a few remarks. 

Firstly, the claimed upper bound 0.6627... is classical. It is related to the fact that, 
as well as in [2, 28], the proof uses the machinery of real-analytic functions. We 
refer the reader to [35, 24] and references therein for general notices. A treatment 
of the argument, as needed in the present paper, is provided in Section A.l. 
Secondly, as it may be seen to the choice of a^, the distances among the planets’ 
semi-axes are not of the same order but grow super-exponentially going towards the 
sun. This resembles a sort of belt arrangement, observed in nature for asteroids. It 
is possible to prove an analogous result, with increasing distances in the opposite 
direction. 

Thirdly, the result in Theorem A (especially, the claimed growth of a^) may be 
regarded as an alternative way of solving the problem of the limiting degeneracy - 
without Birkhoff normal form. 

Acknowledgments I am indebted to Jacques Fejoz, who let me know the work 
by Harrington ]19], without which I would never have thought of this application 
of the P-coordinates. Also, I am deeply grateful to A. Celletti and R. de la Llave, 
for their interest and for encouraging me with precious advices. Thanks finally to 
L. Biasco for his interest. 



CHAPTER 2 


Kepler maps and the Perihelia reduction 


We introduce the Perihelia reduction, or P-map, in the slightly general context of 
Kepler maps. 

Fix a reference frame Go = ) in the Euclidean space E^. We identify 

the three chosen directions k^^'^, k^^'^, k^^^ with the triples of coordinates with respect 
of the system of coordinates established by themselves: 

Definition 2.1. An ellipse (with a focus in the origin and non-vanishing ec¬ 
centricity) is a quadruplet 2 = {a, e, N, P), where a G 1R+ is the semi-major axis, 
e G (0,1) is the eccentricity, N G nS^ is the normal direction and P G N-^ fl S'^ 
is the perihelion direction. 

Definition 2.2 (Kepler maps). Given 2n positive “mass parameters” rui, • • •, 
m„, soil, • • •, DJtn, a set X C K®", we say that a map 

K : K = (Xk,^) G D := X X T” ^ (j/K, xk) G C := K(D) C (R^)” x (R^)" 

where 

e= ,in), iyK,XK) = 

j = l,---n, 

is a Kepler map if there exists an injection 

tk : Xk G X —>■ €k = (<Si^K, • • ■ 7'2 „,k) 
which assigns to any Xk G X an n-plet ((£i,k, • • • , ^u.k) of (co-focal) ellipses 

‘Si.K = (aj.K,ej,K,j = 1, ••• , n 
and K acts in the following way. Letting := x Pj^\ then 

( 2 . 1 ) x^^ = aj^KP^''P^3 ,kQk^ ylP = P^j.kQk^ 

where, if Cj.K) the eccentric anomaly, is the solution of Kepler’s Equation 

(2.2) (yx — ej,K sin 
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then 

(2.3) 

aj,K := aj,K(cosCj,K — e^.K) 


S,K — 




sinC: 


.J,K 




^°j,K — 


'l-e^RsinC: 


.J,K 


iMjil - el^) 


cos (j 


.j,K 


aj,K 1 — ej,K cos Cj,K 

Remark 2.3. The definition implies that 

(i) K is a bijection of the sets D and C; 

(ii) the angular momenta and the energies^ 


aj,K 1 — ej,K cos 


(2.4) 


pO) ._ Jj) y ,,(i) 


jjO) WvkV _ 

^ 2mj ||a;W|| 


do not depend on ij and are given by 


(2.5) 




2 \jyU) 


ttU) _ . 

XIr — - 


2a 


i.K 


(iii) the couples {y ^, ) verify the system of ODEs 


( 2 . 6 ) 


m,- 




'‘l.K 


^U) _ ,,(i) 

j-^K — Uk 




JJ) 


ll 4 ^P 

(iv) Even though canonical maps (with respect to the standard two-form) have 
a pre-eminent role in Hamiltonian Mechanics, Kepler maps are used also 
in different contexts in Astronomy, where being canonical is not required. 
For example, one can consider the Kepler map associated to the “elliptic 
elements” injection 

TEe« : (a, e, F, i, H) —^ £Ee« 

where a = (oi, • • • , a„) are the semi-major axes, e = (ei, • • • , e„) are the 
eccentricities, P = {P^^\ ■ ■ ■ , p(")) are the perihelia, i = {ii, - ■ ■ , in) are 
the inclinations, H = (Hi, • • • , fin) are the nodes’ longitudes. 

The only known examples up to now of canonical Kepler maps are the 
classical Delaunay map De^ (its definition is recalled in the next Definition 
2.14) and the map Dep [27, 7] related to Deprit’s coordinates [12], which 
is recalled in Appendix E. Below, we introduce a new canonical Kepler 
map. 

Definition 2.4 (perihelia reduction, or P-map). We denote as P, and call 
perihelia reduction, or P-map, the Kepler map 

(2.7) P: P = (Xp,£)eDp=ApxT"^(i/,x)eR3"xR3r» 

associated to the bijection 

rp : Xp = (0, X, A, ac) € Xp ^ (€i, • • - , € Ep = rp(Xp) C 


Here, || 7 ;|| := J denotes the usual Euclidean norm of u = (ui,U2,U3) E . 
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defined by means of Definition 2.5 and Proposition 2.6 below. 
Definition 2.5. For a given (€i, • • • , £„) c i?^x-•with 


and masses mi, • • 

•, m„, aHi, • • •, Mn, define 

n 


(2.8) C«:=m,^ 

/fm,a,(l-e2)iV« 

II 

M 

O 

1 < j < n 


i=3 


be the angular momenta associated to (Sj and the partial angular momenta, so 
that 


(2.9) Ce^ ^ = Ce ^ 

i=l 

are the total angular momentum and the angular momentum of the last ellipse, 
respectively. Define the P-nodes 

( 2 . 10 ) 

f X j = 1 

Vj := I Uj := X P^^'> j = l,---,n. 

I X j=2,---,n 


Finally, define 

Ep := {(((£i,--- ,€„) C F;^x---x£;^) : 0 < e^-< 1, Uj 0 0 Vj = l,---,n}, 


and, on this set, the map 



Tp ^ : (di, • • • , £n) G Ep —> Xp G Xp 

= Tp l(Ep) 

where 




Xp = (0, X, A, d, K)Gm.’^ xR^x R!): 

X T" X T" 

with 




0 = (00) • • • ) 0n-l)) I? = (do, ■ 

'• • ,d„_i) 


X= (xo,*** ,Xn-l), f^= (/^O,*** 


— (-^ 1 5***5 -^n) 

defined via the following formulae: 

( 2 . 11 ) 


Qj-i 


Xj-i ■■= 


Z := 


•A:(3) 



a,,(,3){k^^\vi) 




dj_i := < 

1 



pU- 

1) 


{ 0!pU 

-1) (nj_i,Pj) 

G := 

lisg’ 

II 

J 

0 := 

OgO) (Pi.ni) 




Kj-i := 




II 


1 

(ToO) 



j = 1 

2 < j < n 

j = 1 

2 < j < n 


Aj IXHj v^tTij Uj. 
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Proposition 2.6. Let Xp be the subset of ML x R" x R" x T” x T” defined 
by the following inequalities 


>XU + X? - 20? + 2^(x? - 0?)(x?_i - 0?) cos 15. < A. 

(2.12) (xi-i - Xi,i5i) (0, tt) 0 < Xn-i < A„ i = 
and 

(2.13) |0o| < Xo |0i| < min(xi-i,Xi) i = I,-■ ■ ,n - 1. 

The map Tp ^ is a bijection o/Ep onto Ap. The formulae of the inverse map 
Tp : Xp = (0, X, A, §, k) G Dp —>■ £p = (£i,p, • • • , (£n,p) G Ep £j_p = (a^.p, e^^p, 
are as follows. Let ti, • • •, 6„, ii, • • •, i„ G (0, tt) be defined via 

(2.14) 


01—1 . 07 

cost,- = ——, cosij := —I < j <n 


Xj-i 


Xj-I 


{with Qn '■= 0, so that i„ = |) and Ti, • • •, T„, Si, • • •, S„ G SO(3) via 

(2.15) Tj := R3(i5j)Ri(tj) Sj := R3(Kj)Ri(ij), 1 < j < n 
and let 

(2.16) := TiSi • • • T,_iS,_iT, (x,_ifc(3) - x.S.T^+ifc^^)) 
with Xn ■= 0, so that 

(2.17) 

y X?-i + X? - 20? + 2y/(x? - 0?)(x|-i - 0?) cos"^ 




j = I,... ,n- 1 


Xn-l 


j =n. 


Then Cp ^ = Cp ^ o rp and 


a,- p = —— ( — ) 


Cj-.P = 


(2.18) 


\ 


1 - 


lie) 


A? 


iV«= 


0 ) 


lie 


(l)l 




Remark 2.7. 

(i) From Cp-* = orp, (2.4), (2.5) and (2.25), there follows that Cp'* = 

Xp'^ X 2/p^ 

(ii) Pp^ A Np \ Indeed, using the definitions, 

. P^) = xi-ifc^"^ • - T,+iXifc(3) . (fc(3)) 


= Xi-icostj - Xjcosij+i 
= 0 . 


(in) := o rp = = Xi-iTiSi • • • T,_iS,_iT,fc(3). 


We shall prove that 
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Theorem 2.8. The P-map preserves the standard 2-form 

n n 

A dxp^ = (d0i-i A d'di-i + dxi-i A + dAi A d£i). 

j=l i=l 

Remark 2.9. Actually, we shall prove a finer result: the change '■= De£~^o 
P which relates the P-coordinates to the classical Delaunay coordinates (see the 
Definition 2.14) is homogeneous-canonical (compare Lemma 2.16). 

Proof of Proposition 2.6. The formula for aj^p in (2.16) is immediate from 
the definition of Aj. Postponing to below that Cp^ := o 7p has the expression 

in (2.16) (in turn this implies (2.17), the formula for and the one for ej^p in 

(2.18) ), we check that the image set Tp ^(Ep) is included in the domain Ap defined 
by inequalities (2.12), (2.13). From the formula for eyp in (2.18), we have that 
conditions 0 < e^^p < 1 for all j = 1, • • • , u correspond to relations in (2.12). Note 
that the first condition in the second line of (2.12) is equivalent to ej^p 1, as one 
sees rewriting 

(2.19) 

||C«f = - 0? - 

Next, recalling the definitions of 0o, Xo in (2-11), and noticing the relations 

0^. = j = 1, • • • , u - 1, 

we immediately see that conditions ^ 0 ^ ni imply (2.13). We have so checked 
what we wanted. 

Now it remains to check the formula for Cp ^ in (2.16) and the one for in (2.18), 
for any Xp G Ap. To this end, we consider the following chain of vectors 

(2.20) 

^ ^ P(i) ^ ^ ^ p(^') ^ ^ p(") 

11]^ . l^j TLj . 

where ui, ni, • • • , i/„, n„ are the P-nodes in (2.10), given by the skew-product of the 
two consecutive vectors in the chain. 

We associate to this chain of vectors the following chain of frames 

(2.21) 

Gq —^ Fi —^ Gi —^ ■■■ —^ Fj —>■ Gj —>■ Fj_^i —^ ■■■ —G-n, 

where Gq = is the initial prefixed frame and the frames, while F^, 

Gi are frames defined via 

(2.22) F, = (p„ ., S«) G, = (n„ •, pW) j = 1, • • • , n. 

By construction, each frame in the chain has its first axis coinciding with the 
intersection of horizontal plane with the horizontal plane of the previous frame 
(hence, in particular, Vj T and ny T P^£'>). Denote as T^ the rotation matrix 
which describes the change of coordinates from Gj_i to Fy and as Sj the one from 
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Fj to Gj. The matrices Tj, Sj have just the expressions claimed in (2.14), (2.15). 
This follows from the definitions of in (2.11). Then we have the following 

sequence of transformations 


Ti Si 


Go -)> Fi Gi Fj -)> Gj -)> • • • -)> G„ 

connecting Gq to any other frame in the chain. From this, and the definitions of 
the frames (2.22), the formulae for Pp^ in (2.18) and 


=X,_iTiSi---T,_iS,_iT,fc(3) 


(?) (?) (?) 
follow at once. Hence, also the formulae for Cp , which is given by Cp = Sp 

with = 0. 


□ 


For the proof of Theorem 2.8, we shall use three auxiliary maps, that we shall 
denote as P, De£ and De^. The map P is very closely related to P; De^ and De£ 
are well known: in the literature they are often referred to as (two variants of) 
Delaunay maps. 

The map P. Define the set 



with j = 1, • • •, n and 


n 


(2.23) 



Dehne a map 



.n 


■+ 


0 = (00, • ■ • , 0n-i) d = (do, • • • , d„_i) 


X = (xo,-- - ,Xn-l) 


(^0, * * * , ^n—l) 


R — (Ri, • • • , R„) 


r = 


= (ri,--- ,rTi) 


with 
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via the following formulae: 
v(i) . a:(i) 


Rj — 




11^0) II 


Xj-i — 11^0^1 




Kj-i — q;„o) (pj, n_,) 


j = I,-- - ,n 
J = I,-- - ,n 


Sg> • fc(3) 


0,_i - 


jO) 


xii 1) 


1^,-1 = 




Q: x(j-J-) (nj-iI ) 
lIxU-l) II 


j = 1 

J = 2, 


Lemma 2.10. Let Dp be the set o/(0, x, R, i?,r) G M" x 


”xR” x]R"xT"xT"x]R" 
sueh that (0,%, i?, a) satisfies (2.13), and let Tj, Sj and C- '* be the funetions of 
defined in (2.14)-(2.16), with (0,x, replacing (0,x,i?, «)• 

The map P“^ is a hijection from Cp onto the set Dp. Its inverse map 

P: (0,x, R, e Dp -)> (2/p,a;p) e M” X R” 

has the following analytical expression: 


(2.24) 


:=L,TiSi...T,S,fc(3) 


P := X P 

Wp p ' ^ p p 


1 < j < n 


Moreover, the following relation holds: 
(2.25) 


eg) =cg)oP = xg) xyg). 

Proof. With similar arguments as the ones of the proof of Proposition 2.6, 
but replacing, in the diagram (2.20), with sg\ Pp^ with ■pgj]]' nodes 

Pfc, nj; with Pfc, n^, one finds the formula for a;~ ) in (2.24), the formula for 
:= ^ P — 

and hence the formula for 


Sg) := Sg) o P = Xi-iTiSi • • • Tj_iSj_iTj/c(3) 


eg) o p = s^^') - si^'+^) = e^^') 

c p p p 

.(i) ; 


being just the formula for eg' in (2.16), with (0,x,i?, k) replaced by (0,x, d, k). 
With the same argument as in Remark 2.7(ii), we see that xg^ _L egg Finally, the 
formula for yg) is found taking for yg) the unique vector verifying 


0) _ ^p 


R, 


„(i) y ,,0) _ pO) 

iL/p /N Up - '^P 


Lemma 2.11. P preserves the standard Liouville 1-form: 


(2.26) 


+ Xj-idKj-i + RjdTj). 


1=1 


1=1 


□ 
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The proof of Lemma 2.11 uses the following easy lemma: 

Lemma 2.12 ([7]). Let 

X = R3(0)Ri(z)x, 1 /= R3(0)Ri(i)y, C:=xxy, C:=xxy, 

with x,x,y,y . Then, 

y ■ dx = C ■ + C • k^^^di + y ■ dx. 


Proof of Lemma 2.11. We may write 


4'^ = TiSi 


1/4 = 

= TiSi • • • TjS^yW, = TiSi • • • TjSjC(^) 

where 





j(l) ._ 


j = !,■■■ ,n-l 


yb) 

CO 

+ IcL) X /c(3) 

Tj 

(2.27) 

Cb) ■- 

Xj-iS] 



with Xn ■= 0, S„ := id. We also let, for 1 < k < j < n and 1 < i < n — 1, 
ci^'^ = Skifk+iSk+i ■ ■ ■ = TkSk ■ ■ ■ := 

n n 



E d"’. 



Ed”’ 



c**: + s 


(t+ij 

i+1 


m=j 


m=j 


where the product T^+iSfe+i • ■ ■ is to be replaced with the identity when k = j. 
We have the following identities (implied by = Y^k=j 


(2.28) 

n n 

S? = E Cf ^ Sf = ^ Cf = 

k=j k=j 

Applying Lemma 2.12 repeatedly and using (as it follows from (2.27)) 

y^^^ ■ dx^^'^ = Rjdxj, 

we have, for 1 < j < n, 

fc=i 

where, as in the proof of Lemma 2.10, 7), i^ denote the functions Lj, ij in (2.14), 
with 0i, Xi replaced by 0^, xu Note that we have used d!„ = 0, since, by definition, 
in = ^. Taking the sum over j = 1, • • •, n, 

n n n 

1=1 1=1 1=1 

In view of (2.28) and of the definitions in (2.14)-(2.15), we then find (2.26). □ 
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The map De£. The map 

De£: (H,r,R,h,g,'^ ^ ^ ^ 

is defined on the set 


^3n ^ ^3n 


D 


:= |(H,r,R,h,g,r) = (Hi,--- ,R„,hi,--- ,h„, 

gi,--- ,g„,7i,--- ,7„) € X T^” X R!;: : Tj > 0, > 0, 

V j = 1,--- ,n| 
via the following formulae: 


where 


:= R3(hj)Ri(i,)x^, y^Jl := R3(hj)Ri(i,)|/^ 

Dei ^ Dei^ ^Dei ^ ^^^Dei 


~ H • 

ij := cos“^ -J- S (0,7r) 

r 1 


^Dei 


x'~ := rj cosgjfc^^^ + r^ singjfc^^^ 


■= ^ singj)/c(i) + {Rj singj + cosgj)k^‘^\ 


set 


Lemma 2.13 (Delaunay). De£ is a bijection from the domain Dg^ onto the 


CdIS - {(y,^) = (2/('\--- ,2/("UW,--- ,x("))GR3n^K3„ 

nj := X y£ 0, ^ 0 V j = 1, - - - ,n| 

where is as in (2.23). The formulae for the inverse map 


— 1 


De£ : (?/,a:) G Cgr« -t (H,r,R,h,g,T) G D 


are 

(2.29) 


Hj = 


Ti := licg^ll 


hj := OfeO) (fc(i\ Uj) I gj := a^(j) (n^-, x^^)) 


T>el 




Finally, De£ preserves the standard Liouville 1-form 




We omit the proof of Lemma 2.13, which may be found in classical textbooks. 


"1! Mi 
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The map De£. 

Definition 2.14 (Delaunay map). Let 

•^De^ ■ — 5 — (^1,*** : 5 L^ , • • • , Lyj, , • • • , , hi, • • • , hy^, 

gi, • • • ,g„) e : r, > 0, ^ < 1, A, > 0 

^ j 

V j = l,--- ,n| 

and let Eoe^ be the set of n-plets , Gin) where <£j = sat¬ 

isfies 

0 < Cj < 1, rij := X ^ 0, V j = 1, • • • , n. 

Fix positive numbers 9 Jli, • • •, 9 Jl„, rui, • • •, m„. Define 

TDe^ : XDe£ := (H, r, A, h, g) G Xde^ —t ^De^ = (£l,De£, ' ' ’ , ^n.Def) 


with £j,De^ = (aj,De^, 


ArO'). 



= R3(hj)Ri(z,)fc(3) = R3(hj)Ri(*,)R3(g,)fc(i) 

where ij := cos“^ 

We call Delaunay map the map 

(2.30) De£ : Del = (H, F, A, h, g, i) G D^ee ^ {yuee, x^ee) e x K^" 
which is defined on the domain 

Duei ■= ^'Del X T" 

as the Kepler map associated to roe^ via the following lemma (the proof of which 
may be found in classical textbooks). 

Lemma 2.15 (Delaunay), roe^ is a bijection of X^ei onto Eds^- Its inverse 

map 

'^De£ ■ — (Gll.De^i * ' ' ; Gln.De^) ^ Ei3e£ ^ Xi3e£ G Xj^e£ 

is defined by equations 

rH, = c«.fc(3) r r, = iic«ii 

(2.31) < < Aj = mjytMjOj, 

[ hj := afc(3) [ gj ■= a^u) (n^,P(^)) 

where is as in (2.9). Furthermore, Del preserves the standard 2-form 

n n 

^ ^^dL = ^ {fH^j A dhj + dVj A dgj + dAj A dlj). 
j=i i=i 

Now we are ready to complete the proof of Theorem 2.8. 

Proof of Theorem 2.8. Let 

D* := {P = (0,x,A,i?,r,£) G Dp : P(P) G Coe^}. 
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It is enough to prove Theorem 2.8 on Dp, since indeed the P-map is regular on 
Dp = Dp. On Dp, we consider the map 

:= Der'oP: 

P = (0; X) £ Dp Del = (H, P, A, h, g, £) G Djjg^ := ipDeei^p) C Doe^- 


gives the Delaunay coordinates at left hand side in (2.30) in terms of the 
P-coordinates at left hand side of (2.7) in the subset Dp of Dp the P-image of 
which lies in the De£-image of Dper- Clearly, leaves the (A,£) unvaried. More 
precisely, decouples into two disjoint maps: the identity on the (A,£), and a 
4n-diniensional map 

'■ (0iXj^;^) G Dp —>• (H,r,h,g) G Dpg^ = </'Der(Dp) 

on the remaining coordinates, which turns out to be a bijection of the sets Dp and 
^Dee- map and the sets Dp and D^ei do not depend on (A,£). 

Indeed, the explicit expressions of Dp in terms of P = (0, y. A, i?, k, i); or of 

T>jyet in terms of Del = (H, P, A, h, g, £) involve only the Cp \ the P^^^, 

that do not depend on (A,£): (2.31) (where one has to replace C with P), (2.15) 
and (2.18). 

In view of the previous consideration and of Lemma 2.15, Theorem 2.8 is implied 
by 

Lemma 2.16. The map preserves that standard 1-form: 

n n 

(Hjdhj + Tjdgj) = Y + Xi-id^i-i)- 

t=i i=i 


Proof. We look at the analogous map 


b^~ 

Der 


(0,X,^,«) G D| ^ (H,r,h,g) G D^^ = cPYi^l) ^ D 




'Def 


The analytical expression of this map is identical to the one of This follows 

from the fact that depends on the coordinates (0,x, i?,k) only via C~ ^ and 


p 


exactly as depends on (0,x, d, k) only via Cp^ and lip'', that C 


r(i) 


■•(i) 


,0) 


and — fjy- have exactly the same expressions of Cp and pbO, apart for replacing 

Ibp II 

(0,X5'd, k) with {Q,x,'d,K). Compare (2.29) (where one has to replace C^^ with 
C~^), (2.31) (where one has to replace Cg^ with Cp^), (2.15), (2.18), (2.24) and 

(2.25). But Lemmata 2.11 and 2.13 imply that (p~^ preserves that standard 1-form: 

n n 

Y + fjdgj) = Y {^3-iddj-i + 

3=1 3 = 1 


Then (p^^^ does. 


□ 
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2.1. The P-map vs rotations and reflections 

Now we discuss how the P-map behaves in presence of symmetries in the Hamil¬ 
tonian due to rotations or reflections. 

Let H = H(y, x) be the Hamiltonian governing the motion of n particles, where such 

particles are expressed in the canonical coordinates 

Assume that H is left unvaried by rotations and reflections. Namely, if 

0R.S : ^ {Ry^^\Sx'^^'^), j = l,---,n 

where R, S are real a 3 x 3 matrices, then rotation invariance is 

H o = H V R : RR* = id 

while reflection invariance is 


i 

( cn 

0 

0 \ 1 

f n 

0 


H o = H for some S^- = 

0 

0-2 

0 s. = 

0 

T2 

0 

\ 

V 0 

0 

CT3 / 

V 0 

0 

7-3 / 


CTi, Ti = ±1. 

Rotation invariance is associated to the conservation, through the motion, of the 
total angular momentum is (2.23). Reflection invariance is not associated to 
integrals. 

The Hamiltonian Hhei in (3.1) is rotation and reflection invariant, and reflection 
invariance holds with any choice of cr, r. 

Let 

Hp :=HoP. 

The fact that ^ is preserved along the motions of H implies that the coordinates 

©0 = Z, -do = C) kq = 0 

do not appear in Hp. Indeed, Z and C are integrals, while g is conjugated to 
G = ||Sp^||, which is an integral for Hp. Thus, the number of degrees of freedom is 
naturally reduced by two units, once one regards G as a prefixed external parameter. 
Namely, for any fixed xo = G, Hk may be regarded as a function of the 2(3n — 1) 
dimensional coordinates 

P := (0,X, A,d,«;,£) 

which does not depend on kq. Here, 

0 = (01, • • • , 0„_i), d = (-di, • • • , -dra-i). 

An analogue property is also shared with the action-angle coordinates (d>, P, A, 
Ip, 7 , t) described in [27, 7], and related to a set of coordinates discovered by A. 
Deprit [12] in the 80s (compare also [36] or the Appendix E). 

The main novelty introduced by the P-coordinates (that does not hold for the 
coordinates of [7]) is how P behaves relative to reflections. 

We denote as 

R.^ := (('S_^(2) .S^(2) 0-^^) = (1,-1,1) 

the reflection of the second coordinate both for the y*--^Ps and the a;*-'^Ps and we let 
S"(0,X, A,-d,K,^) := (-0,x,A,-d,K,£). 
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Proposition 2.17. 

(2.32) R^oP = PoS". 

Therefore, ifil = ll{y,x) satisfies 

H o R- = H 

then Hp := H o P satisfies 

Hp oS" = Hp. 

Hence, any of the points 

00 = • • • = 0 „_i = 0 , {tto,-■ ■ = {ko,-■ ■ ,kn-i)Tr mod 27rZ" 

is an equilibrium point for Hp, for any {x, A, k, €). 

Proof. Defining := T^Sj, := Tjk^^\ we write the vectors Pp^ and 
Sp^ (compare Eq. (2.18) and Remark 2.7(iii)) as 

= r(i) ... = Xj-iR^^^ • • • 

The explicit expressions of R^-l^ and s^^'^ are 

(7) 

R)]^ = cos Kj-i cos ttj-i — sin Kj-i cos lj sind_,_i 

(i) 

R 21 = cos Kj -1 sin '&j-i + sin Kj -1 cos ij cosd_,_ 1 

ii) 

Rjt = sinitj-isinij 

R^'^^ = — cosij sinKj_i cosi?j_i + sindj_i(—cosij cosij cosKj_i +sintjsinij) 

R 22 ^ = — cosij sinKj_i sindj_i — cos dj_i(—cosij cos cos +sintjsinij) 

( 7 ) 

1^32 = cos ij COS Kj -1 sin tj + sin ij cos tj 

Rj'^^ = sin ij sin Kj_i cos i?j-i + sin dj_i (sin ij cos ij cos Kj_i + sin ij cos ij) 

R 23 ^ = sinij sinKj_i sindj_i — cos'dj_i(sinij cosij cosKj_i + sinij cosij) 

R 33 = — sin ij cos Kj_i sin ij + cos ij cos ij 

s^^'^ = sinij sindj_i 

82 ^ = — sinij cosdj_i 

U) 

S 3 ' = COSij. 

_ (7) (7) 

Then S lets Pp and Sp respectively, into 

(Pj(j))- : = R-and ■=-Rf ^. 

Therefore, (with := 0) and x P^^'> are 

l|Cp II 

transformed, respectively, into 

(C«)- := -R^-C^, (q(,^))- := RfQ<^\ 

On the other hand, Oj^p and Cj.p are left unvaried by S“. In view of Definition 2.2 
and Definition 2.4, the thesis (2.32) follows. □ 




CHAPTER 3 


The P-map and the planetary problem 


After the reduction of the invariance by translations, a Hamiltonian governing the 
motions of n planets with masses ^mi , • • •, /iTO„ interacting among themselves and 
with a star with mass mo can be taken to be the “heliocentric” one 


(3.1) Hhei:= ^ 


l<i<7 


V 2mi 


vciiMi 

II.T(b|| 


M E 


y(d . yio) 


mo ||xb) — II 


where (y, x) = • • ■ , • • • , are “Cartesian coordinates” taking val¬ 

ues on the “collision-less” phase space x R^” \ A, where 

A = |a; = • • • , G R^ X • • • X R^ : 0 ^ ^ V 1 < i < j < n| 

endowed with the standard 2- form 


n 3 

:= dy A dx := ^ ^ A 

i=i j=i 


and with 
(3.2) 


dJli = mo + ymi 


mom^ 

=- 

mo -I- ymi 


being the so-called “reduced masses”. 

In the following Section 3.1 we describe a general property of Kepler maps, in 
relation to their application to the Hamiltonian Hhei- Then (in Section 3.2) we 
shall specialize to the case of the P-map. 


3.1. A general property of Kepler maps 

For a general Kepler map K, we denote 


Hk(K) := HheioK = -^ 


^ 2aj_K(XK) 


m/k(K), 


where 


/k(K):= ^ 

l< 2 <j<n 

and y^, a;|^^ are as in Definition 2.2. 
We denote 


{i) {j) 

vk • yk 

mo 


mimj 


ii4’-4"i 


(3.3) 


/k(Xk):= 


(27r)'^ 


fK0^K,i)d£, 
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SO that 


/k= E /k. 

l<z<j<n 


fK= Y. fK 

l<2<j<n 


— 
Jk •— 


(i) (j) 

Vk ■ Vk 

mo 


mimj 


„(*) 


,(i) I 


fij _ ^ 

Jk ■— 


(27r)’^ 


fildii 


■ din 


For a general Kepler map, one always has, as a consequeirce of (2.6), 


1 


/ Ti^'^dL = ^ 


1 


27T , 


27T 


(j) 

jp . _ /Ti(i) yO') _ _ 

aij — -r Vy — 


mm 


2a 


i.K 


(3.4) 




1 


,. 0 ) 


2^ ' di,-0, 


U) 


where we have denoted as 


rpU) ._ WVK^ ip 

K ■ 2m, 


IT kc 


T/(i)_ 

K - 


the kinetic, potential part of in (2.4), respectively. 

Consider the average /k(Xk) in (3.3). Dne to the fact that ylp has zero-average, 
one has that only the Newtoniair part coirtributes to /k(Xk): 


/k = - Y 


miTTij 


diidij 


(27r)2 Jj2 II 


T2 


l<2<j<n 

We now consider any of the contributions to this sum 


/•b _ _ 
Jk — 


mimj 


dijdii 


(27r)2 

and expand any such terms 


/t 2 ||a;^^ - 


1 < i < j < n 


—^ ^(0) ^-(1) --(2) 

fi'3 _ fV _i_ f^3 _|_ 


where 


—^{h) 

fij 

Jk 


fJ — f- 
Jk — JK 

mimj 


+ /k 

1 d^ 


/T2 


h\ de^ 


(d 

r* ^ 


IJ) I 


e=0 


diidij 


is proportional to ^(f^)^- Then the formulae in (3.4) imply that the two first 
terms of this expansion are given by 


-^(0) _ m^mj 

JK ~ 7 ■ 

aj,K 




Namely, whatever is the Kepler map that is used, the first term that depends on 
the secular coordiirates Xk is the double average of the second order term 


/k (Xk)--7^Zw4^- 


(27r)2 


Now we specialize to the case of the P-map. 
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3.2. The case of the P-map 

We denote as 

(3.5) Hp(Xp,f) = hO„,t(A) + /i/p(Xp,£) 
where 


Xp := (0,x, 


(A) 


(3.6) 


h° 

^fnst 


■ 1 ’ 
j=i J 


is the Hamiltonian (3.1) expressed in P-coordinates. 

Using the definitions, it not difficult to see that 

Lemma 3.1. /p , /p depend, respectively, only on the coordinates 

Xp '■= (©i, • ■ • ;0jA(n —l)j Xi—Ij'"' jXjA{n—l)j A^, Aj, di, • • ■ l)j — l) 

W:= 

with a Ab denoting the minimum of a and b. 


2 , 2+1 
P ’ 


with 


Accordingly to the previous lemma, the “nearest-neighbor” terms f 
i = 1, ■ ■ ■, n — 1, depend only on 

(3.7) 

(©i, ©z+1, Xi —1: Xi^ Ai+1: Ai, A^^i, tti, Ki^ 1 < I < TT _ 2 


x; 


2 , 2+1 


(©72—1: Xti—2: X 72 —1: Ayj_l, A^, ttn—l, 


i = n — 1. 


2 , 2+1 


( 2 ) 


we have a special rule. Indeed, for any Kepler 

^J2) 


However, for the functions /p 

map K, the “exterior” angular momentum is an integral for . This 

r(2) 


readily implies that any is integrable, for having four degrees of freedom 

and four independent, commuting integrals (||Cp^ -|- (Cp^ -I- Cp'*’^^) • 

I|Ck~''^^|| and itself). This fact has been firstly noticed, in the three-body 

case (* = 1, n = 2), by R. Harrington [19] who, using the Jacobi reduction of the 
nodes Jac, where the coordinates are named 


G„ 




A,: 


f- 


i = l, 2 


(with Gi = ||C*'®^||, gi related to the perihelia directions, and G := ||G||, C = 

G*^^^ -I- GG) appearing as an external parameter), noticed that depends only 

on (G, Gi, G 2 , 7 i, Ai, A 2 ). 

r(2) 


Let us now inspect how the integrability of /p*^^ is exhibited in terms of the 

7 ( 2 ) 


P-map. Since ||Gp"^|| = Xti-I: one has that /p does not depend on Kn-i, 

and hence, by (3.7), depends only on 


X n — l,n 
13 


■— 1 , Xn —27 X.n—lj 
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This fact, for n > 3, is no longer true for z = 1, • • •, n — 2, because in that case 
Xi l|Cp“^^^|| (indeed, x^ = However, since, for (0i+i,z?i+i) = (0,7r), 

||Cp''’^^|| reduces to 


lie 


d+i)| 


(0i+i,i5i+i) = (O,'n-) 

one has that the functions 


= Xi - X*+l 


=( 2 ) 


^2,2 + 1 _ r2,2+l 

/p /p 


( 2 ) 


i = 1, 


i = 1, 


n — 2, 


n — 2 


(0i+i,'i3i+i) = (O, 7 r) 

do not depend on Ki and hence, by (3.7) depend only on 


Xp ■— (0i, Xz—1; X*5 XZ+1; -^2+1; ^ ‘ ^ 

In the following lemma we provide their explicit expressions. 


Lemma 3.2. The function fp 
=(2) 


and , for n > 3 and 1 < i < n — 2, the 


functions /p*^^ have the following expressions 
(3.8) 


_ aLi A3 r5 

/p — .35 ry 

4a3 xi-i L2 

3 40Li-xL 


(30^_l Xra-l) 


I2 ^ {xl-2 + xl-1 - 20L1 + 2y^ (xLi - ^l-i)ixl-2 - ©Li) cosz?„_i 


, 3 ixl-l - Ql-i){xi-2 - ©Ll) . 2 Q 

+ --T 2 -sm dn -1 

and 


Jp 


f(2) 


= mirrii+i 


A3 


^302-x2) 


4af+i X?(Xi-Xi+i)^ ^2 

3 40f-xf 


(3.9) 


2 A? 

3 
2 


(x*-i + Xi - 20,^ + 2y^(x* - 0D(xLi - 0?) COST?* 

- 2 

sin di 


3 (xf-ef)(xti- 0 f )_2 


A? 


Lemma 3.2 is proved in Appendix B. Here, we limit to the following. 
Remark 3.3. 


(i) The formula in (3.9) holds also for complex values of the coordinates, 
provided that 

arg(Xi - Xi+i) e (-f > f] mod 27 r. 

(ii) The importance of the formulae in (3.8) and (3.9), which is the main fea¬ 
ture of the P-map, is that, exploiting the equilibrium for (0^, idi) = (0, tt), 

_(2) l'2^ 

the integration of /p~^’” and of /p*^^ can be performed explicitly, 
switching to a suitable associated convergent Birkhoff series, as Lemma 3.4 


below states. Direct integrations of f 


n—l,n 


(2) 


for example, starting with 
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Hamiltonian computed in [19], appear technically much more involved 
and, up to now, are not known. 

Lemma 3.4. It is possible to find complex domains ‘Si with non-empty real part 
and a canonical, real-analytic change of coordinates 


(3.10) 

where 




y* := 


Ji ■= 


(3.11)x, 
such that 


{X*n-2,X*n-l,K-l,K) 

(k* k* f* 

iXn—27 X**—1; A^) 

(Xi-i.Xi,X*+i, A*, A*+i) 

(^n— 2 5 ^n—1 j ^n— 17 ^n) 

7 ^i7 "^*+1) 

1 - 1 .^ j,n-l 


i = n — 1 

i = l,'--,n — 2&n>3 
i = n — 1 

i = 1, - ■ ■ ,n — 2&n>3 
i = n — 1 

i = l,---,n — 2&n>3 
7 = n — 1 

i = 1, - ■ ■ ,n — 2&n>3 


/, 


^ J.'hl 


(3.12) 
depends only on 
Y* 


h* *= 

^^SCC • 




:( 2 ) 


i = n — 1 


z = l,---,n — 2 & n>3 


/ Pn-l+gn-l A* A» * » \ 

V 2 ") Xn—2’> Xn—l) 


i = n — 1 


( P^i +‘>i A* A* V* V* V* 1 
i 2 I ’ ^'-*+1’ Xi-ii Xj j Xi+il 


2 ’ ’ ^'-*+ 1 : Ai-iJ A* j Ai+ 1 / i — l,---,n 2&n>3. 

T/ie transformation may be chosen so as to verify 

y*=y*, (©*, x^-^i) =F^{p^,q^,y*) 


(3.13) 

if 


^ini(-K’-9*,yr,x*) = (-0,,-l?i,yi,Xi) 




Lemma 3.4 is proved in the following Section 5.2.1. 









CHAPTER 4 


Global Kolmogorov tori in the planetary problem 


In this section we show how the P-map can be used to prove Theorem A. We 
defer to the next Section 5 more technical parts. 

4.1. A domain of holomorphy 

A typical practice, in order to use perturbation theory techniques, is to extend 
Hamiltonians governing dynamical systems to the complex field, and then to study 
their holomorphy properties. 

In this section we aim to discuss a domain of holomorphy for the perturbing function 
/p in (3.5), regarded as a function of complex coordinates. We shall choose it of 
the following form 


Dp := Te+_.,j+ X (Xe x T^) x (Ag x T^), 
where, for given positive numbers 

0+, d+, G±, A±, s 

with i = 1, • • •, n, j = 1, • • •, n — 1, 

(4.1) 

Te+.,,+ :={(0,d) = (0i,--- ,d„_i)GC"-i xTJl-i : 

|^,-7r|<d+, |0,|<0+, y j = !,■■■ ,n-l] 


Xe :={x = (xo,--- ,Xn-i) e C" : < |Xj-i - Xj I < G+, | Im (x^-i - Xj)l < 

V j = I,-- - ,n| 

Ae:={A = (Ai,-.- ,A„)gC": A-<|A,|<A+, |ImA,|<0, 

V j = I,-- - ,n| 

Ts :=T + i[-s,s] 

with Xn ■= 0. 

The domain Dp will be determined as the intersection of the “collision-less” set, 
where, as functions of complex variables, the mutual distances of the planets 

dj,P ■■= Ibp^ - II 
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are far away from zero, with the holomorphy domain of P, where, again as as 
functions of complex variables, the absolute values |ej,p| of eccentricities in (2.18) 
are bounded away from 0 and 1, those of the inclinations |tj|, |ij| in (2.14) are away 
from 0 and, finally, Kepler equation (2.2) provides a holomorphic solution. 

The latter issue is not a peculiarity of this problem, since it naturally arises in 
the context of the two-body problem’s equations. In the early XX century, T. 
Levi Civita [24] studied the holomorphy of the solution of Kepler’s Equation with 
respect to the eccentricity. The holomorphy with respect to the mean anomaly has 
been investigated, using similar arguments as in [24], in [4[. Here, we address the 
problem of determining the holomorphy with respect to both the arguments. 

Proposition 4.1. Let e = 0.6627... be the solution of 

(4.2) 0<p<l & / -. = !■ 

1 -I- v'T + p^ 

Then for any 0 < e < e, one can find a positive number £ depending on e such that, 
for any e = ei -|- ie 2 G C, with jej < e, the complex Kepler equation 

^ — e sin ^ = £ 

has a unique solution C{£,e) which turns out to be real-analytic for £ £ T^. 


The following result completes the study of the holomorphy of /p. 


Proposition 4.2. Let e be as in Proposition f.l. For any given e^, Ci, with 


0 < < Ci < e i = 1, ■ ■ ■ ,n 

it is possible to find positive numbers 

A„ B„ C. >C„ d" , sG(0,l), a £(0,1) 

such that, if the following inequalities are satisfied 
(4.3) 

C,A+ < G- < G+ < C,A-; 

n-l 

^G 


max 




et 

Gn 


dj < min 




< 5 


then the eccentricities inclinations u, \i and the mutual distances di^p verify 


(4.4) 


e* < |ez,p| - 


,p| < Ci, max||costi|, |cosij||<cr, |di,p| — 


Proposition 4.1 and Proposition 4.2 are proved in Appendix A.l and A.2, re¬ 
spectively. We shall use them in the form below. We remark that the super¬ 
exponential decay of the semi-major axes ratio will be used only in Section 5.2 
below. 
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Corollary 4.3 (choice of parameters). Fix e, < Ci, c G (0,1), and let C, < 
c: < c; < C,, D, := min{Av'(S )2 - (CJ 2 , B^(Q )2 - (C*)^}, D := mini<,<„_i g 
a < ^. Define, for i = 1, - ■ ■ ,n and j = 1, - ■ ■ ,n — 1, 

(4.5) 

G+:=C:A-, G-:=CrA+, 0+:= sG", 

'^n 

6i 

where af is as in (*). Then fp is real-analytic in the domain Sp. 


4.2. A normal form for the planetary problem 

Definition 4.4 ([ 2 ]). Given m, ■ ■ ■, Pm G '■= + ■ ■ ■ + Vm, let 

£o D £i D £2 D • • • D = {0} 

be a decreasing sequence of sub-lattices of IF defined by 
(4.6) 

2,0-.= IF, 2i := [k = {ki,-■ ■ ,kyn) Gl'', kj G : ki = ■ ■ ■ = ki = 0} 

with i = 1, ■ ■ ■, n. Next, given 7 , 71 , • • •, 7 m, t G K+, we define the set 
of the (71 • • • 7 m; T)-diophantine numbers via the following formulae: 


D 


,,, .= jueK-: |u-t| >vt e i;,_i \ft, |t|, < A'} 


.K, I 


- n d; 


2=1 


i',K 

7l-"7m;T I I -^7l-"7m;T- 

kgn 


with kj € , 


(4.7) \u;-k\^ 


i=i 



forany fc = {ki,- 

• • , km) G Z‘^\{0}, 

71 

\k\- 

if fci 7 ^ 0 ; 


72 

\k\- 

if fci = 0 , ^2 7 ^ 0 ; 



> 


. \kn 


if ki= ■■■ = km-1 = 0 , 


, kjn 0 . 


Remark 4.5. The choice m = 1, 71 := 7 gives the usual Diophantine set Dl( 
The TO = 2-case, with 71 = 0(1) and 72 = 0(p), where p is the strength 

of the planetary masses has been considered in [2] for the proof of the Fundamental 
Theorem, mentioned in the introduction. 

The following result in proven in the next Section 5. It is unavoidably detailed. 

Proposition 4.6. Let m^, be as in (3.2) and := J — 

2, • • • ,n. There exists a number c, depending only on n, mo,-- - ,mn, a^, Cp ej, 
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and a number 0 < c < 1, depending only on n sueh that, for any fixed positive 
numbers 7 < 1 < K , a > 0 verifying 

c 


(4.8) 
and 

(4.9) 


K < 


a‘ 


.3/2 


1 r ,a+,5i?2f+2 K2{f+i)^ 

- ^<1 

c l. 7 ^ 7 ^ J 

there exist natural numbers i^i, - ■ ■ , i' 2 n-i, with Uj = 3n — 2, open sets B* C 
B'^.,X* C X, positive real numbers > ■ ■ ■ > 72 n-iei, • • • , ^n-i, • • • , f„_i, ri, • • • , r„, 

a domain 

®n := B^ X X- X Ap X x Tf, 
a sub-domain of the form 

D* := B*^xXtxArX x 


measD^ > (l — —) meas:Dr 


verifying 

(4.10) 

a real-analytic transformation 

K- (P, 9, X, A, K, £) G S*-)> 2)p 

which conjugates Hp to 

Hn(p,g,X, A,K,£) := Hp 0(()„ = hfa^t,5ec(P,<?,X,A) + /i/eyp(p,g,X, A,K7 
where /epp(p, g, x, A, k, £) is independent of kq, and the following holds. 


1 . The function hfasi,s^ciPjQ,Xi-^) *■5 a sum 

hfasi,sec(p, q, X) A) = hfast(A) + fj. 9, X> A) 

where, if 


y^ ■= 


pf+qf pl_i+ql_i 


2 ’ ’ 2 
then h^dst h^ec given by 


Xi—lf ‘**5 X . n — 1 -, 


h,...(A) = - E ^ - f E 


A? 


h«c(p, q, X, A) = ^ h;ec(yi) 


where the functions h* have an analytic extension on Sn o,nd verify 


, _ 1 («^,)^ 


(a„_j+i)3 


< mec(yi)l < 


C (a„_j+i)3 



-cK 


Oi 


2 . T/ie function /epp satisfies 
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3. //C is yi deprived ofxo, the frequency-map 


C t wjo5t,scc(C) •— c(C) 


is a diffeomorphism of x Xi x A~) and, moreover, it satisfies (4.7), with 

m = 2n — 1, T = f > 2, and 


1 
2 


Uj := 3 


j = I,-" ,n 
j = 3, n = 2 
j = n + l, n>3 
n-\-2<j<2n — 2, n>4 
j = 2n — 1, n > 3 


(4.11) 


9a, h 


fa0t,scc 




ojj := 


d hfast.sec 

(— --,Xn-2,Xn-l) 


d 

( 


^2n —j4"92„_j 
2 


,X2,,-3-i) 


lfast,scc 


9p2 + ,2 hjust 
2 


J = I,-" ,n 
j = 3, n = 2 

j = n + 1, n > 3 

n-|-2<_)<2n — 2, n>4 

J = 2n — 1, n > 3 


' ±1 
a- e, 

M(a+_„)^ 7 

(Oj + l-ra)^ %-r!, 


1 < J < n 


n-|-l<_)<2n—1 


4. T/ie mentioned constants are 
£j := c rj 


with T > 2. 


K^+i' 


ri := c 6j 


4.3. A “multi-scale” KAM Theorem and proof of Theorem A 

In this section we state a “multi-scale” KAM Theorem and next we show how 
this theorem applies to the Hamiltonian Hn so as to obtain the proof of Theorem 
A. 

Theorem 4.7 (Multi-scale KAM Theorem). Let m, ui, • • • , Um G N, u := ui-f 
• • ■+Vm > ^, r* > u, 7 i > • • • > 7 m > 0, 0 < 4s < s < 1, pi, • • • • • • ,ei> 
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0, Bi,--- ,B( C Dj := eR: {x,y) e Bj} C M, B := x • • • x C 

D ■= DiX ■■■ X C C. A:= DpXCr. Let 

H(p, q, I, V') = Hp, 9. 1) + f{P, 9: 1, V”) 

be real-analytic on B^^ x Cr x T|^^, where \i{p,q,I) depends on (j),q) only via 


J{Piq) ■= ( 


p\ + ql 


pl + qj 


)■ 


Assume that ojq := is a diffeomorphism of A with non singular Hessian 

matrix U := d^j(^p g) /)h and let Uk denote the + • • • + Vm) x v submatrix of U, 

i.e., the matrix with entries {Uk)ij = Uij, for vi-\ -hi^fc-i +1 < i < 1 < J < 

where 2 < k < m. Let 

M>sup||J7||, Mfc > sup ||C/fc||, M > sup ||C/“^||, E>\\f\\p^s+s 


Mfc > sup ||Tfc|| if U '^ = 
A 



1< k< 


Define 


fEMlL\ 

\~ir) 


where log^ a := max{l, logo} 


s 

Pk ■— 3]\/[^^T»+l ’ ^ iPli ■ ■ ■ ! Pmi Pi, ■ ■ ■ ) Pi, Xi, ■ ■ ■ , r,y—c} 


L := max 
EL 


{m, Mr\ ••• , M-i} 


E := 


32 ■ 


Then one can find two numbers c^, > Cu depending only on v such that, if the 
perturbation f is so small that the following “KAM condition” holds 

CvE < 1 , 

for any w G O* := uJo{D) fl one can find a unique real-analytic embed¬ 

ding 

(pa,- -d = 0,^) €T’' -a + 

G Rea X X Re5^ 

where r := c„Ep such that := is a real-analytic u-dimensional H- 

invariant torus, on which the H-flow is analytically conjugated to d ^ t) -\-ujt. 
Furthermore, the map (diuj) -A fitjfd) is Lipschitz and one-to-one and the invari¬ 
ant set K := (^ T(^ satisfies the following measure estimate 


meas ^Re (a)xT"\K^ < Cj^^meas xT")+ meas (Re {Dr)\D)xT”‘^, 
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where denotes theuJo-pre-imageofY>j^^...^j^^r^ inD. Finally, onT^ x 

f 2 », the following uniform estimates hold 

kfc(-;w) - 

|u(-; w)| < c,yEs 


where Vk denotes the projection ofv= {v,v) G x • • • x over R^*", Vk '■= 
and J°(w ,I^(uj)) G D is the uiQ-pre-image of uj G O*. 


|WfcP 

2 


Theorem 4.7 generalizes [6, Proposition 3] in two respects. The first generalization 
concerns the consideration of m > 2 scales (in [ 6 ] only the case m = 2 was treated). 
The second consists of taking H depending also on the rectangular variables {p, q) G 
Such generalizations can be easily obtained, and hence will be not discussed 

here. 

Proof of Theorem A. Let 

^ ■.= c\/a(\oga~^y'^^, .A=ilog — 

c a 

where c is as in (4.10) and c will be fixed later. We aim to apply Theorem 4.7 to 
the Hamiltonian Hn of Proposition 4.6, with these choices of 7 and K. To this end, 
we take 


cia^. e] 


M, = <! 


/r(a^ 


+ 12 


. ci(a^+i)30 


^ < 3 


n + l<j< 2 n —1 




C 3 a. 


Pj = 


C 5 


C 5 




KT,+i 
7^J —n 


E = 


1 1 ( 02 )^ 


1 < j < n 
n + l<j< 2 n —1 


ce 7^ (ai )^ 


g-cA^2(r.+l) 


r iCr 1 ^*1(4)^ 

L = M =- ' 

C2 f 


K = 


1 




C 4 °^V 7 ^(a-) 


« 0i7 

p := - r* > 3n — 2 

Kr, + 1 


where K := max{A, A}. The number ^ (a-^s bounded by ^ for a suffi¬ 

ciently large N depending only on n. Hence, if c < and a < Cg, we have A < 1 
and the theorem is proved. □ 















CHAPTER 5 


Proofs 


In this section we provide the proof of Proposition 4.6. This is divided in two 
steps: normalization of fast angles and of secular coordinates. 


5.1. Normalization of fast angles 

— — (fc) 

Let /p , fp as in Lemma 3.1, and let 

—(>2) —^ --(0) 

(5.1) . 

Proposition 5.1. There exist two small numbers't, Ci, where'c depends only 
on n, while Ci depends only on n, mi, • • • ,m„, such that, if the inequality in (4.8) 
and 

(5.2) 

hold, one can find a real-analytic and symplectic transformation 

4>fasi ■ (0) X) ^; •^) £ ®fast := T^e+,?i?+ ^ ^ X T~, X T~, —£)p 

which conjugates Hp to 

Hfa5t.erp(0, X) A, d, K, £) := Hpoc/ifosi = hfasi(A)+p /fa5t(0, X) A, d, n)+p /fa<ii.erp(0, X) A, -d, k, £) 


where hfa^t is as in Proposition j.d, and 

n—1 

(5.4) /fast := /fast: /fost.cfp 

i=l 

Here, 


E/f-i ,epP' 


1. The ‘fast frequency-map” 

to^fast ■= ^h^dst 


is a dijfeomorphism of A with non-vanishing Jacobian matrix on A^g and, moreover. 


Wfasi G D 


^ asi 
7fast,T 


V Ag A, 


with 

7fast := (7i> 
and Vi, 7 i as in (4.11); 


7«) Pfast := (Pi, • • • , Pn) 
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2. the functions ffast ejp depend on kq; the given by 

(5.5) 

-^(>2) ~ 

/fast =/fa<ii(ti,y*,Xi) =/* (ti,yi,Xi) + i = l, 

with 




j=i+l 


‘ ‘ ; Qn—1; ^2; ‘ ‘ ; ^n—l); 

yt ■” (Xi—1) ■ ■ ■ ) Xra—1) -^2) ■ ■ ■ ) .^n) 

Xi .— (^2? * * * 5 ^n—l)- 


/n particular, rfo noi depend on l\, - ■ ■ 


3. finally, /expjast satisfy the following bounds 


(5.6) ^ 


y ^.i_i 


-^2+1 


II /fast,cyp II® fast — 


1 e 


—'cKs 


<^1 “i+i 


Let Lq, • • •, L„ be defined as £i in (4.6), with v = m = n and pi = • • • = = 1. 

Lemma 5.2. If K verifies the inequality in (4.8), then one can find a number 
C 3 , depending only on mo, • • •, mn, such that 
C 3 


|^k,fast(.^) ‘ ^1 ^ 


-)3/2 


y k G Lj_i \Lj, |fc| < K, V A £ Ag, V j — 1, 


Proof. For A e Ag, Wk.fast.j := verifies < Iwk.fasi.jl < 


(a-)3/2- 


In the case j = n, we find [wk.fasi ■ k\ = |wk,fasi.nfcn| > ’ since A:„ ^ 0. Let 

then j ^ n. For k £ Lj_i \ L^-, kj 0, so, inequality (4.8), with C 2 < 
and (4.5) imply 

and hence 

l^k.fast * ^1 


miiij .y/gq - 
maxa \/9Jl 


K < 


min, 


_ _ liiiii I I I 

max,- ^ af ) 


in 

:n-l V 


a^\ 3/2 


= I y^Wk ,fasl.ifc 2 | > inf |wk,fast,j| “ AT max sup jwk,fasl,i 

« ^ An f <' r<i n. \ 


> 


j<2<ra As 

f'—J 

V^, - maxi>, i/IOti v^,- 

— AT—^— > -^— . 


(a+)3/2 (a-, i)3/2 - 2(a+)3/2 


□ 


Proof of Proposition 5.1. The proof proceeds by recursion, in n steps. 
We describe the h^^ step of this recursion, with h = 1, • • •, n. We start with a 
Hamiltonian of the form 


(5.7) 


Ifli-l ~ + p, fh-l 
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where is as in (3.6), and a domain 

®ii-i = T0+(h-i)_.,9+(h-i) X Xff(h-i) X Ag(h-i) X X T”(h_i). 

When h = 1, we take Hq := Hp, 0^°^ := 0+, := i9+, := 6», := s, 

fo ■= fp Eind we decompose 

n —1 n 

fo ■■= fo ■■= ?h with /* := ^ fpf 

i—1 


We observe that /g depends on the coordinates 

*** 5 0n —15 Xi—lj *** iXn—Ij Aj, 

1^2; * * * ; 1; ‘ ‘ —15 * * * ; ^n- 

For n > 3 and 2 < h < n — 1, we assume, inductively, that fh-i is a sum 
(5-8) fh-l = fh-l + /erp,?i-l = fh-l A /epp,Ii-l) 

l< 2 <n l< 2 <n 


where, in turn. 


with /^_i depending only on the coordinates 

: ©n—15 ''' ; -^n 

' ' * : 1 : ^25 ‘ ' 5 ^n—15 ^iWhi ‘ * 5 

and /^_i, /^_i 5 /cyp,/i-i verifying the following bounds and identities 


fLi = 

\\fLl\h.-^ < 

(5.9) ll/cVp,i2-ib.-i < 

Here Hp^ denotes the projection over the module L/j. In any case, /i = 1, or 
2</i<n—1, we focus on the Hamiltonian 

n—1 ___ 

(5.10) H,,_1 = h^„^t + fifh-l = h^ast + ^^Y1 fh-l- 

i=l 

Our purpose is to apply Proposition D.l to this Hamiltonian, in the case that the 
abstract system (D.l) does not depend on the coordinates {p,q). To this end, we 
take the coordinates 

I := A, (/?:=£, 77:=(0,x), ^:=(fd,K), 


ai 
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the functions fi in (D.3) to be the /^_i, and 


N = n — 1, V = rrii := 2i 

{^1} ' ' ■ : tVi) 

(^15 ‘ ) ■ (^nj * * * 5 ^max{n—) 

(^1; * * * ; Vfrii) *—(071—1;*** ;0n—7;X77—1;*** ;X77—7—l) 

(Cl; * * * ; ^rrii) •— (*^n—1; * * * ; ^n—1; * * * ; ^n—i) 

^7 ■ — (-^77 ;*** ;-^l;071 —1;*** ;077 —7;X77—1;*** ;X77 —7—1;*^77—1;*** — ^n— 1; * * * ; — 


The non-resonance assumption (D.2) for uj = Wk.fast = ^Ahk^fast; with 
37 = La-1; 3 = U737 = La-i, £ = La K = K 

is ensured by Lemma 5.2, with 

A n{h-l) 


a = 


(a+)3/2 


A = A, r = e\^ 


Now we have to check condition (D.7). In the case 2 < h < n — 1 the inductive 
assumptions (5.9) and assumption (5.2) imply 

ll/Cb._7 < WlLilh.-i + WfLilh.-^ < (i + )^)lliCb.-2 

(5.11) < * * * < (1 + G:i.A-ici)'*-i7^||x,„ < =: E,. 


An analogous bound holds also for h 
evaluated as 

Ci = e(l -b 2 ze )/2 di = min {0 


= 1. The numbers Ci and di in (D. 6 ) may be 

(/ 7 -l)^(A-l )^0 + (/ 7 -l)^ + (/ 7 -l)| ^ C2e[^-^\ 


From these bounds it is immediate to see that inequality (D.7) is implied by (5.2), 
provided Ci < 2 “^|(|)"'“^C 2 /(£ 4 C„). Then Proposition D.l applies. Its thesis 

implies that Ha-i in (5.10) can be conjugated to a suitable H)) = Lk,fasl + M/a; 
where verifies equalities and inequalities in (5.8)-(5.9) with h replaced by h -|- 
1 and £i a-i, ^ 2 ,h-i replaced by suitable £* a; ^2 h- Then, applying the same 
transformation to Ha-i in (5.7), we shall conjugate Ha-i to Ha = hkjast + p/a; 
where /a satisfies the same equalities and inequalities as /^, with suitable £i,a > 

cLa. ^2,h > e:s,A* 

After we have performed n steps, we let Sjost := '^n, Hjast^ejp := H 77 , := 

^asi /fast~/p; /fosi.ejp *” /exp,n; /fast := X]7=l /*> /fast := X^fcl // /fast.cfp : = 

Efc/ /fai;t.efp - ^ith /* := E"=7-Pi fp' Therefore, 


Hfast — + /^(/fast + /epp.fost) — 


(0) 


l< 2 <j<n 


fp /fast H“ /exp,fost) 


reduces to (5.3) and the formulae given below, using (5.1). 

It remains to check the bound on the left in (5.6) (the one on the right follows by 
construction). This follows by telescopic arguments. Indeed, 
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ll/fa^lb. 


< 


/pb„ = ll/A - = l|nL„/* - nupy. 

n 

^||nL„/^-nL„r^iCbn 


h=l 


< 


< 


ElinL„/^-nL„nur^Cll 

h^l 

n 

/i^l 

n 

EllH-nL.T^iCb. 


/i=l 


“ ^ a7 aj 


■^i+l 


Here, we have used the second bound in (5.9), (5.11), that does not depend on 
ii, ■■■, £n, and, finally, nL„ = nL„I> = nL„nL,., for all 1 < h < n. □ 


5.2. Secular normalizations 

Consider the following truncation 

Hfast(0,xA,^,K) := fifast(A) +/r/f„,5i(0,x,A,d,K) 

of the Hamiltonian Hfasi,efp in (5.3). The purpose of this section is to describe 
an iterative scheme which, after (n — 1 ) steps, conjugates Hjost to a close-to be 
integrable system, with an arbitrarily small remainder. 

Let us firstly establish the following notation. 

• Given a Taylor-Fourier expansion of the form 

9{p,q,K)= E 

(a,b)eN2”^ 

we denote as 

^p,q,K9 ^ ^ 90,a,a 


m 


9a,b,k 


p-\q 

V2 


P + iq 

V2i 


{p,q,K) e H2™i(0) X 


Proposition 5.3. There exists number Ch, depending only on n, mo, ■ ■ ■ , mn, 
such that, for any h = 1, ■ ■ ■, n — 1 and any K, 7 > 0 such that (4.9) hold with 
c replaced by Ch, one finds open sets 

G*CG,:= 

verifying 

(5.12) meas [B* x G*) > Tl — meas(H^^. x Gj) 

such that, defining 


G7 




j = n-h, 


n — 1 
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(5.13) 


rph _ 

— 


< 


\^j -TT\<Ch-^, |0j|<c?,G+ 

V J = I,-- - l| 


n>3, 1 < h < n — 2 


otherwise 






X 


IfeS.chf •— {x — (xo, • • • , Xn-i) : Xi-I - Xi G iG*)-c^0., Xj-1 - Xj G (GjOchP^- 
V z = 1, • • • ,n - - 1, j = n - h,- ■ ■ ,n, Xn := o} 


:= T-^, X B;^ X X x x T? , 

a real-analytic transformation 


‘^5ec,h • ®scc ^ ®fast; 

may be found, which conjugates ffast io a new function 

/scc.Zi •“ /fast ® *^’sec,?i 

enjoying the following properties. 


1. Denoting by where 



t(0 


= (0f\ •• 

c)(^) .,9(^) 

5 '^n-h 



= (pW,<7W) = 

(Pn-h^ ■■ 

• 0^"*^ • 

’ Pn—1') ^n—h'> 

■■ > 9i-\) 


y(0 

= (x(0,A(")) = 

=(xr, •• 

, An-1’ ^'■1 ’ 

'• , aW) 

(5.14) 

x(0 

= (k('*),^W) = 

••• 

(^) ... 


coordinates 

on t/zen 

zs co-variant with the symmetry: 



4>sec,Zi(-t^'*\-Z 


(0) = (-t(O),y(O), 

xG)) if 




yW^xW) 

= (t(°),y(°),x(0)) 


and hence, 

fstcp 

1 is even around 





= (0,fc7r), 

z(^) = 

: 0 k€{o, ly 

1—h—l 
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2. Defining 


:= 


A%-i) i<n-h-l 

0 otherwise 

/(pfV±(9fV (h) 


'■W 
Yi = 




2 ’ 
A(^) . . . 


5 Xz—15 ' ‘ 1 \n—l'i 


i > n — h 


ipj^i^f+ieio^ (h) ^ 

2 5 Xi—i’t 1 Xn—li 

otherwise 


-(h) 
x; = 


j • • • 5 n^n-h-^) n> A k, l<h<n — 3Sz l<i<n — h — 2 


otherwise 


(5.15) 


and y := yi, x := xi, fsec,h has the form 

= K^cAfnlh) + 

(5.16) + /efp,5ec,/i(t^'*\z('*\y('*\x('*)) 


KAfnlh) = Yl ^«c(yf^) 


(5.17) /„o.m./.(t("\y("\x(0)= ^ fAmMiA\yr\^r) 


3. the functions /norm h decomposed as 


Kccin i = HcciJi i^HcAn ) 


-I- 

1^0. lo; /normA'^^^ ’ ~ /normA'^^^ )^Jnovm.hy^i ’Xi ’ ) 


( 5 . 19 ) 


fnotm,h — Y^ n?i(/] 
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and “ Lemma 3.4- The functions /norm,ft, /crp,scc,ft (5.16) 

may be bounded as 


\Kc.h\ < -max{AiK(^) 
C/j i. Oi 


an\3/2 1 


'‘i+1 


\P hi < 

1*^ norm,All — 




• max 




i+l 


(5.20) |/ejp,acc.ft| < 

Cft (a„ y 


{atf gf+i (a+)^ 1 

7 (o+pi)^’ ^i+i (a/+i)^'’ 

{at? -I 

7 (a-+i)3/ 


4. Defining 

( {Pn-hf + (^i-ft)^ i.Pn\? + {(ln-1? , (ft) 

V 2 , ■ ■ ■ , 2 ’ Xi-i 

so i/iai 

= aW) 

/or any ■■■, h-^n \ the map 

:=9c('»hoec,ft(C^"\AW) 

is a diffeomorphism of Dr x with non-vanishing Jacobian matrix. The set 
D^ X Xi consists of the subset of Dr x X^ such that wjast.scc G > where, if 

Vj, are as in (4.11), 

ft'scc •” (Pn+1, '** , P2n—l) 7scc (7n+l,*** , 72n—l)- 



We shall give the complete details of the proof of Proposition 5.3 along the following 
sections 5.2.1-5.2.4. In this section we just provide main ideas. 

Scheme of Proof. The proof is by recursion. The hf'^ step of this recursion 
starts with 


/sec,ft—1 — hjec/j—i -f /norm,ft—1 T /cfp,scc,ft—1, 

where, for h = 1 

(5.21) hjcc 0 = 0, /efp,sec,0 = 0, /scc,0 /norm,0 ■— /fast, 


while, for n > 3 and h = 2, ■ ■ ■ , n—1, we assume, inductively, that hsec_ft_i, /sec,ft-i 
and /epp,sec,ft-i satisfy the theses of Proposition 5.3, with h replaced by {h — 1). 

The transformation conjugating /sec,ft -1 to /sec, ft will be constructed as a 

product o '/'norm of an “integrating” and a “normalizing” transforma¬ 

tion. 


Due to the bound on /cfp,sec,ft-i, it is enough to focus on the truncation 
(5.22) 

n—h 

_T, if _ u I \ ^ fi 

Jstc,h — 1 ■— ^sec,/i —1 “T /norm,/!—1 — ^^sec./i—1 i / ^ Jnozm,h—l\^i ) 










5.2. SECULAR NORMALIZATIONS 


45 


of fscc,h-i- We split 

n—h—1 

f , _ fu-h (Ah-1) Ah-1) Ah-l). ^ „ (Jh-l) Ah-1) Ah-1). 

Jnovm,h—l J novm,h — l\^n—h ^ n—h '>^n—h /”*” / ^ Jnovm,h—l\^i "i J i ^ ) 

and we distinguish two cases. 

Case n>3, /^ = 2, •••jn—1. By the inductive assumption (see (5.15) with h 
replaced by (h — 1)), the function depends only on 


^n-h — \^n-h > ^n-h ) yn-h 


therefore is integrable. In Section 5.2.2, we shall construct a canonical, real-analytic 
change of coordinates 


:n-h 

nnl 


Tlh T)^-^ 

"^inl ^ "^scc 

t-W (h) (h) (h) 




(5.23) DL := . x B?^ ^ x x X*’)r^- x Ai,r x T? , x T? 

^ Chtf ChSn-h ChS C/iC',C/ir ^h' ChS Ch 

such that 

(5.24) 




norm.Ii-l " V^int — if *,n-h) 

depends only on where is defined analogously to y-^^ in (5.15). Here, 

tl"^ := (0l"\ #)) 
yl") := 

Jh) ( (h) n(h)\ 

are defined analogously to (5.14). 

We shall construct in such a way that it involves only the coordinates 


t(h-l) ^Q(h-l)^ ^(h-l)^ 
^(h-l) (^p(h-l)' q(h-l)^ 

yih-1) _ \x^h-l)^ A('‘-b) 

x(^-l) := £ih-i)^ 


. ((h) (h) (h) . , Ah-1) (h-l) (h-1) 

'^inl ■ y*,n-h^ ^*,n-h) ^ Hn-/t ’ ^n-/t+l’ J n-h ’ ^n-h ) 


with 


„(0 

^*,71—h 

{h) 

y*,n-h 
,Sh) 


= ) 
\P^,n—h'> ’ 1 ’ ^*,n—h'> ’ ^n—lj 

— A/ A ('ll! 

VX*^rt-Ii-l> ’ X*,n-1) rt-Ii) ) ■‘'■n J 


h) _ (Ah) 

,n—h ’ \^*,n—h—l^ 


(h) Ah) 

’ ^*,n—1’ ^*,n—h'> 


(5.25) 


'^n-h ■— \'^n-h > ^n-h J 

Ah-A — (Ah-A ... Ah-A Ah-A ... 

^n-h+1 \fn-h+li i Pn-1 i 'Jn-h+1^ i 'Jn-1 J 

Ah-A (A^~^'> ... A/ ... A 

yn—h ■ \Xn—h — l^ > Xn—l^ ^In—h ’ ’ ) 

Ah-A .^ ( (h-A ... ^('*-1) fih-A ... f(h-A) 

^n — h ■ \^n—h—li ’ ^n—1 J ^n—h ’ i ) 
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and has the form 


^n—h int /i’ * / 


(5.26) CT" : 


^n—h '^int v^*,n—/i’ ^^,n — h'> / 


'',3 


(/i-l) (h) 

yn-?» =y*,u-h 


. ^n-h^ =^i%h + V’[’^Upi%h^ <li%h^ y*'"^) 

with f|^], odd, (p[^l even in 


{pi^l-h+i? + {(li^l-h+if (pi^i? + i<li^i? 


fi^'’ ■■= ( 

dh) _ ((h) (h)\ _ (h) j_-(h) 


’ ^,n—h/ 


(5.27) gi:^):=pl>igi:^ 

with j = n — h + 1, ■■■, n — 1, ioi n > 3, h > 2 and as in (5.25). 

In particular, observe that <('’^7^ enjoys the following properties: 

• it is co-variant with the symmetry: if 

yl"\ xl^)) = (t^'*-!), y(''-i), x^'*-!)), 

then 

yi'^\ xl^)) = (-t('^-i), -z('‘-i), y(^-i), x('^-i)); 

. leaves the “actions” 

yi^^ = 

unvaried, where is as in (5.27), and 
-(/i-l) _ (iPn-h+l) +(gn-/i+l) (Pl ) + Wl ) ,.(?i-l)\ 

y ■” I 2 ’ " ' 2 ’ 

is defined analogously; 

• leaves the averages with respect to the x-coordinates unvaried. Namely, 
for any real-analytic function g on 


Applying to f^cc.h-i in (5.22), we obtain 


n—h— 1 


/scc,int,/i—1 •— /sec,/i —1 ' 


_ -U _i_ _L \ ^ 

iJ-sec,/7 —1 “T ^^scc ^ Jnovm,int,h—l\^4‘,i ^ y"> 


2=1 


2—/l—1 


= V h* V r /.(^) 

/ ^ ^^scc,h\y ^ J norm,int,/i— 1\^*,2 5J*,2 5 -^*,2 / 


i=n—h 


2=1 


with 


(5.28) hsec,/i := hsec.?i-i + hi 


i—h 


fnotm,mt,h—l * /norm,/ 2 —1 


n—h 
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and (as it follows from (5.15) with h— \ replacing h and (5.26)) depends 

only on the arguments 


(5.29) 



Ah) . 


, PfA) Ah) _ _ _ Jh) \ 

’ ^*,n— /i — 1’ ^*,2 ’ ’ ^*,n— /2 —1/ 

~ih) 

y*J '■ 

— (pA) 

12 I (Ah) ^2 (Ah) n2 , (Ah) ^2 

^{h) \P*,n—l) ' W*,n—17 

2 ’ ’ 2 


A.*,^— 1 1 

A/ 

^{h) , 


, ■■■ , n>Akl<h-l<n-3 


1 0 

otherwise 


The next step will be to retain the dependence on 

and, for /i < n — 1, to to eliminate from /scc,ini,ii-i the dependence upon the angle 




to an exponential remainder. Namely, we look for another canonical, 
real-analytic change of coordinates 


(5.30) 


l^n-h 


^scc ^ ^int 


(t 


(h) 


ZW, ^(h)^ 


(Ah) Ah) (h) (h). 


so as to conjugate /sec,int,/i-i to a new Hamiltonian 
(5.31) 


^—h— 1 


p ^ _ p in—h _ 1.^ \ '' pi ^ ^ (^)’\ _j_ £ 

Jsec,h ■— Jsec,ini,h-1 O 'Pnotm ~ + 2_^ J notm,hAi jJi i j + /efp.scc, 


i=l 


where /*„„„ and /ejp, 5 cc,ii satisfy (5.18)-(5.20). We choose as the subset of 
where the map 


(p';!'lA^+(dhA (h) h = 2, ■■■ ,n-2 kn>4 

2 ,Xn-h~l 

^ ^Hcc.ra-i h = n — I 
2 

does not verify resonances up to order K, and next we apply a suitable normal 
form theory (Proposition D.l). We shall choose (t>n^zm such a way that 

• it is co-variant with the symmetry: if 


y^'^), x('')) = (tl^\zl^\ yi’^\ 




then 


(5.32) Co™ (-t(^), -z^'^), y^'^), x^'^) ) = {-e\ -zf), y 


(h) Ah) Ah)\ 


• leaves the “actions” 


.ih) 


(h) 
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unvaried, where 

nC'*) 'iS , (Jh) 


^n—h ’V o ’ o’ 


... 

^n—h'> ’ An—1’ 


K\ •••, AW); 


„(h) . ^(pi.„-?.+i)^ + (9i.„-/.+i)^ (p 14 T + (CiT ^(h) ^{h) 

y*,n—h ■ (, 9 )■■■ 9 ) X*,n—?i> '"' 5X*,n —1) 


2 


(5.33) ••• , 

• verifies 


.x)'*) . , , (5 ° Cotm) = , 5 ) o ' 


(5.34) nr.) 

+ + l ' .. ' ^Ti-h,+ l’"^n-h + l 

The thesis of Proposition 5.3 at rank h follows, with 

'i—h 


in—h 
^norm ■ 


— h 


/sec,/i ■— /sec,/! /cjp,sec,/i —1 ^ ^sec ’ /cpp,scc,/! ■— /cpp,sec,/! “1“ /epp,scc,/!—1 ^scc 

Case h = 1. The proof of this case uses similar ideas as the proof of the case 
2 < h < n — 1 {oi n > 3. However, due to subtle differences between the two cases 
(compare, e.g., the inductive assumption on in (5.15) for h > 2 with Eq. 

(5.35) ; the definition of h"“^, for h > 2 in (5.24), with the definition of h"“^, 
in (5.37) and (5.40)), for sake of precision, we briefly discuss also this case. 

Let /scc.o be as in (5.21). In view of (5.4) and (5.5), we can split 

n—2 

-^ 

(5.35) 




£ _ £ri—L,n- j_ rn—1 

Jscc.O — Jp ' JP 


/;A + E/h 


fast 


where 


- -(>3) -^(>2) --(2) 

xn—1 ._ rn—1 _ rn—1 

J P ' «/ P «/ P 


—(2) 


and the summand appears only when n > 3. As for /p ’ , by Lemmata 3.4 

(see also Lemma 5.4), we find a domain lB„_i (defined in Eq. (5.41) below), a 
real-analytic and canonical transformation 

(5.36) 

. cd 1 d 1 (C yT _1 

^int • y*,n-H ^ 45„_1 —>■ ^ X)„_l .— 

of the form (5.26), with h = I (but neglecting the coordinates zf \ such that 


(5.37) 

depends only on 


/p"-i^'^c;' = h?c7\yi!Li) 


(5-38) yi^Li = 
In (5.36), we have let 


f (1) (1) ,,(1) ,,(1)^ 

[ -^- 2 -’ X,,n- 2 : x;.n-l, A). ;,_p Ai,4J- 


(Ji) o(i) ) 

v(l) I 

y*,n-l •“ \A*,n- 2 i A*,n-l5 
^( 1 ) (M) .( 1 ) \ 


r p(0) /0(O) n{0) ^ 

^n-1 ■— ^n-l) 

AO) / (0) (0) ^(0) ^(0)N 

Jn-l \Xn-2^ An-1’ ) 

AO) (AO) (0) X 

k -^n-l •“ V^n-2’ ^n-l)- 
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We let 


' t(0) := 

y(0) _ (;^(0),A(0)) 

^ x(0) := (k(0),£(0)) 


:= (0l^\ 4 ^^) 
7.4 := {p4, q4) 
:= 

„(i) ._ (M) /lit 

-A.* .- 7^* ) 


analogously to (5.14), with h = 0,1, and then we regard the map in (5.36) as a map 
CT^ : (tl'\ yi'\ xl^)) (t(°), y(°), xW) 

on the set 

-Oint ■— I ( 1 * : I y* I j ■ ( 2 *,n-lJ y*,n-lJ ^*, 71 - 1 ) ^ 45 „_iJ 

where is defined on the extra-coordinates via the identity. 3),^^ has the form 
in (5.23), with h = 1. Applying this extension to /sec.o in (5.35) we obtain 

n— 1 

/«c,int,0 := /scc.O o = h”7^(yi^Ll) + /notm.ini,o(ti^],yM,X^^]) 


where 


rn—l , / p 

•/norm,ii^t,0 ' V] 


n— 1 ^—I 

+ /fast ) ° ^int : 


/, 


norm,int,0 ■ f fa6t ^int 


and, as a consequence of (5.5) and of (5.26), with h = 1, q depends only 

on the arguments 

^(1) ._ ("0(1) ... 0(1) .a(I) ... ^(1) ) 

^*,2 • ’ ^*,n—2’ ^*,2 ’ 2’/ 

'(1) .= ('Ji) „(i) y(i) ... y ••• A^^^) 

■y *,2 ’ ^*,n—A*,2—’ A.*,71 — 


i(l) ('^(1) ... ^(1) \ 

,2 • V^*,2’ ’ ^*,n—1/' 


( 1 ) 


Cri.int.o is a function of 


Note, in particular, that 
(5.39) 

(tii!,n—li y*,n— 1 ; X.i;^,.j_l) = (p.i;y^_xj 17 X*, 7 i —25 X*,n—1; 15 


In view of the fact that h"ec^ depends on the actions in (5.38), we aim to eliminate 
from /sec,ini, 0 the dependence on the following angles 

K*.i if n = 2 

^*,n—2: ^*,71—1 if Tl ^ 3 

and to retain the dependence on {p4n-i^ 9*^1- 1 ) ouly via 
we choose a domain 2)5|,e ®int as in (5.13) where the frequency 


(pi(Li)-+(e-i)\ Then 


d ( 1 ) .O I / (1) yO 

p*,n-ir+(K/n-ir ^ (1) 

2 :X*,n-l 


n—1 

•sec 


<^scc,l • — 


^ (1) ..2 I / (h) n2 

2 —2’^Hc.n — 


n— 1 
sec 


n = 2 


n > 3 
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is non-resonant up to the order K and on this domain we construct a real-analytic 
transformation (5.30) which conjugates fsec.i to a Hamiltonian 


/sec,! ■— /scc,int,0 O 


in-1 

^norm 


n—1 

: i:. 

2=1 




(1) ,Hi) ,Hi) 


) + / 


cj:p,scc,l 


Now, since (as it follows from (5.39)), f^otm i i® actually a function of only, 
this step is proved, with 


(5.40) 


Osec IJn-ll ■— Osec lyn-ll + /norm.liyrt-ll- 


□ 


5.2.1. Construction of The following lemma completes Lemma 3.4. 

In particular, it provides the transformation in (5.37). 


Lemma 5.4. Let i = 1, • • •, n — 1. Let A, X, 9 in (4.1) be chosen in such a 
way that 


inf|g|>0, sup|arg 5 |<^ 

Sp Sp 4 

(5.41) V5e|xj-i, Xi, Xi-i+Xi, 5 xj_iA 2 - (xi_i - x,)2(4xi_i - x»)|- 


Then, the domains 05^ in (3.10), the functions h'^,. and the transformations 
can be taken as follows 


= 


^int 


Bl X Al, X xgi X 

Bl X Al, X x|. X T| 
Pz 


i = n — 1 

i = 1, - ■ ■ ,n — 2 & n>3 


(5.42) h*,, = A, 


©i = j: + Mp^,Qi,y^) 
di-TT = Piqt A gi(pr,qi,y*) 

yi = y* 

^ Xi=x* + pi{p^,qi,y*) 

E, + + 0(p„ ft)® 












5.2. SECULAR NORMALIZATIONS 


51 


where x denote the projection of the set Xg x Ag over the coordinates yi in 
(3.11), 9 9/2, s := s/2, fi, pi are 0{pi,q^)^, odd in {pi,qi), ipi is 0{pi,qiY, and 


e,. = y 


A := 




- (x^-l - - Xi) 

+Xt) 


A, := m^TOi+i 


4«i+i 


E, := 

n, := 

Ti := 




2(x* - Xi+i)^ 
_3^f+i_ 

X*Af(Xi - Xi+i)^ 


5-3 


(Xi-I - Xif 

A? 


(5xi-iAf - (xi-i - x*)A4xi-i - Xi))(Xi-i + X^) 
9 (Xi-i - Xi)^(3Xi-i - X*)(5Xi-i + Xi) 


xliXi - Xi+i)^ )■ 16 xl-lX^^^/3t 

.o^ 3 2x/_i+9x?_iXi+ 2x*-ixf+ xf 3x*-iX?,, , ^a4l 

(5.43) -g--_ 


with Xn = 0, Ci depending at most on the ratios a/'/a~, the masses mi, • • •, mn 
and, as usual, 'y/z denoting the principal determination of the m*** root of a complex 
number z. 


Proof. Since the formula for /p coincides with the one for /p taking 
Xn = 0, we shall only work on the terms /p*~'’^’s. 

Let Xi be as in (3.11), and let 

(5.44) D, : (0„ A) e S A^^ x x, e 

where Tq+ is the projection of Te+_, 9 + over the coordinates ( 0 ^, 1 ?^), while m^ 
is 4 or 5, accordingly to (5.42). We shall obtain the transformation in (3.10) 
as a product Ant = ° where Atog ^b\x described below. 

A Taylor expansion of /p*^^ around (0i,di) = (0,7r) gives 


(5.45) 


pi,i+l 

Jp 


Ei + rii- 




(i^i —tt) 


+ R-i 


where A^, E^, Pi, Xti are as in (5.43). Note that Pi, Xli are well defined under 
the assumption (5.41). The expansion in (5.45) shows that (0i,i?i) = (0,7r) is an 

elliptic equilibrium point for /p*^^. The remainder is given by 


R,: = F 


3 402- 
2 A? 




(x?-xti)^ 


+ 2J(xf - 0f)(xLi - 0?)(1 + cos A) + 


(^/W^+^/xll-of)^ 

i(x?-0?)(xLi-0f) 


A? 


sim tli 
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where the symbol F on the left means that only terms of the fourth order in (0i, -di — 
tt) have to be included. The lower order expansion of Ri is 

R^ = n.i0f + T2,i(dj - 7r)^0j + + O(0i, i9i - nf 

with 


n,* 

:= n(yi) 

3(xi_i - Xi)2(3xj_i - x»)(5xi-i + Xi) 

8x?-iX*A2 


:= 7-2 (yi) 

3(2xf_i + 9xf_iXi + 2xi-ixf + xf) 

7-2,* 

4x*-iA2 

'Fa,* 

:= 'Tsiyi) 

Xi-ixh. , ^ 

:= g^2 (4Xi-i+Xi)- 


We introduce the generating function 


c ~ Q ^ 

*hoia0,i(Pi,yi,di,Xj) = -=-hyjXj. 


A 


It generates the canonical transformation 
Pi 


Kiag ■ ~ 


■di-TT = yi=yi, x^ = x^ + 


^y,|3^iy^) 

A (A) 


PiQi 


which transforms into 


(5.46) 

with 


/oiag.i — /] 


*>»+! o rA* — A ■ 

° ^Tiias ~ 




A := /3(A), Ai := Ai(yi), E, ;= C'(A), := ^{fi), 

R* := Rj o = Ti^^pj + f 2 ^^pfqf + Ts^iqf + 0{pi,qif 

n(A) 


Tii := 


f2,^ ■= T2(A), A,i := T3(A)A 


To compute the domain of (/'jjjjg, we use the following inequalities, which readily 
follow from the definitions: 

Ci „+ < lAI < ~ 


and 


|A.A(A)|/ 1 


< 


A(A) CiA’ 

We then see that, choosing a suitable Ci < q, and the domain 

: \{pi, qi)\ < A = A\/A a g Ag X Xg. Xj e T™) 

4 & 

inequalities^ (5.44) are verified, as desired. Now we look for another canonical 
transformation 

Abir: (Pi,ft,y*,X*) -t {pi, ft,A,A) (y* = A) 


^Compare (4.5). 
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defined in a analogous domain 

: M,qt)\ < Ci = c*y* e x 

with c* =: Ci < til2^ such that 


G T” 


fiiag,i ° (Ihn ~ ^sec 


satisfies the thesis of the lemma. We aim to apply Theorem D.4, with 

h = E, + a^^^, / = R„ e = 2c*v^, e = c*v^. 


We have to check that inequalities (D.20) are satisfied. We can take a and e as 
it follows from the following inequalities, which, in turn, are easily implied by the 

definitions 



inf |9h| 
<3* 

= 

inf|H.|>^^^=:a 

sup Ri 
<8* 

< 

sup|R*| < , maxsup{ , (d* Tr)^} 

£>* s* (Gn)^ 


< 

(G)"(G+)2 

-=: c 

Ci 


Here, we have used that, for |(p^g*)| < 2c*v^, (0*,-d*) := ^ipt,qt) veri¬ 

fies 


i0:i 

1^:- 


M 

m 


< 2c* 


Gt 

'i's/^ 



ti 


\qm\< 


2c*v^v^ 

Gn 



We then have that condition (D.20) holds, provided one takes 



From (5.46), one easily computes that the fourth order term of h|(.j. corresponds to 
be as in (5.42), with 


3 * ,1 * ,3 , 



^li ■■= 


Finally, properties (3.13) easily follow from the construction. 


□ 


5.2.2. Construction of (l>\nii " ' •> — 3)- We have to solve (5.24), 

assuming that Proposition 5.3 holds, up to rank h — 1. Accordingly to (5.18), 
(5.19) and letting 


^n—h+l 

^inl 


in—h-\-l 




O (j) 


n—1 

'int 
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we may split 

rn—h _ 

«^norm,/i—1 


^ - r(>2) _ 

E TT f ^ I rn — h 

j—n—h-\-l 

_(2) _ _(>3) _ 

n,.-io + Uh-l(f^-h,n-h+l - ^ 

_ ” -(> 2 ) — ' 

E n,’ 

j—n—h+2 

=(2) 


^n-h+1 


l\ _|_ rn—h 
/ ”*” / szc,h—l 


{^) 

pTl — h.TL —/l-j-l I TT / p'. 

Jp +11/1-1 

_(>3) 

n ( pTl — h.TL —/l-j-l 

h-l[jp 


£71 —h 
5tc,h — l 


( 2 ) 


~n—h,n—h+l ^ 


^int 


■)+ ^ n,_i(/’ 


—r^(>2) 


j—n—h-\-2 


int 


pn—h,n—h-\-l^ — ^^ _ pn—h,n—h-\-l^—^^ pn—h,n—h-\-l^‘^^ 

Jp Jp ~ Jp 


pn—h,n—h-\-l 

J] 


.( 2 ) 


P 

=( 2 ) 


-( 2 ) 

pn—h,n—h-\-l^ ' 

Jp 


pn—h,n—h-\-l 

' J'P 


:C2) 


=( 2 ) =( 2 ) 
jn-h,n-h+i Lemma 3.4. Note that we have used that jg 

unvaried by Let ^n-h, 4\nJ Lemmata 3.4, with the symbols 

-h, -dn-h), Yn-h, ^n-h of that lemma corresponding to 


^int 

); Yn—hi ^n—h 
-{h-i 


^n—h ■ \^n—h ’ ^n—h ) 

.>-1) \ 

yn—h ' V^n—/i—1’ Yn—h ’ A-n—/i+l’ ^^n—h ’ ^^7i—h-\-l) 

^n-h ‘^n-h ’ ^n-/i+l’ *^n-h ’ ^n-h-\-l) 


Defining 


—*(h) ( *(h) 

^TL-h ■= [PtL-H^ 



—*(h) ( *(h) *ih) 

Yn-h ■= [Xn-h-V Xn-h^ 


Y*(^) ^*(^) 

^n-h V^n-/i-l’ ^n-h^ 


X 


*{h) 
n—h-\-l ’ 




*{h) 

n—h-\-l'> 


A *(^) 

^^n-h- 

f*W 

^n-h^ 


K*i.h) N 

.^{h) \ 

^n—h+l) ' 


^*{h) ._ ^0*(/i)^ 
^*{h) 

y*(/l) ._ ^ 

^*{h) ._ 


h+1 


in an alagous way as in (5.14), we regard as a map on the set 
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extended via the identity on the extra-coordinates. We then have that trans¬ 
forms into 


£71—h _ _ £71—h in—h _ y^n—h , £n—h 

Jstc,h-1 Jscc,h-1 ^ ^int “ "t ^scc 


where 


£71—h 

^sec 


(5.47) 


-- ( 2 ) _ 

Tj j fn-h,n-h+l ^ i^ri-h\ 

O'l’int j 

_(>3) _ ^ _(>2) _ 

I TT ( — h,7l —/l+l — h\ I \ -TT / £71 — h,j /FvTl — h\ 

+Ilh-l{fp °^ini)+ 2^ ^h-l[fp O^int ) 

j—7i—h-\-2 


_|_ £71—h 

5ec,h — l 


O (j) 


71—h 

int • 


Here, we have used o that Tih-i and commute and 

observe that has the form of in (5.23), with ih replaced by a suitable 
of the same form. The function f^ec^ satisfies the following two properties: 

• It depends only on 


{P 


^{h) 

7i—h'> 


q 


*{h) 

71— 



where y*(^l is defined analogously to (5.27); 

• is uniformly bounded by the right hand side of the first inequality in (5.20) 
(this follows from the definition in (5.47)); 

• is even for 


yPn-h^'in-h) ^ yPn-h^'in-h)- 


Proceeding in a similar way as we did for the construction of in the proof of 
Lemma 5.4, we may apply Theorem D.4, with 


h = h^e7^ / = f?e7\ (P,Q) 

(p! n'\ - v-v*^^^ 

,w ) — ^n-h^ y —y„_^. 


\Pn-h’ 


q 


*(h) 

n—h 


) 


X = X 


*{h) 

71—h' 


with defined analogously to (5.25) and defined 

analogously to in (5.27). We then find another domain as in (5.23) and 

another real-analytic transformation 

such that 


Hni,h 


£71—h ._ £71—h 

«/scc,/i —1 ■ «^SCC,/l—1 


/n—_ £71—h 

Jsec,h—1 


in—h 

hni 


7—h 

<,int 


h 


71—h 

sec 


as desired, depends only on y, in (5.15), and hence (5.24) is satisfied. That 
may be also chosen of a form analogous to (5.26), with \ 

replaced by p*^^\ z*^0, x*(0 also easily follows from the 

properties bove. Therefore the composition 


iTi-h ^ r. 




has again the form in (5.26) and satisfies (5.24), as wanted. 


□ 
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5.2.3. Construction of •••, ^ 3). In this section we aim 

to determine, for n > 3 and l</i<n — 2, a transformation (t>na-!m solving (5.30)- 
(5.31), assuming the Proposition 5.3 holds up to rank {h — 1) and that has 
been constructed. 

We switch from the coordinates defined implicitly via the right hand 

side of (5.30) to the auxiliary coordinates 


G 


(h) 

auj: 


_ /p(^) 


,G 


{h) . 

auj:,n ) 


„(h) _ , ih) 
g(3uj — iBauf.li 


(h) 

I gauf,i 


defined via the linear transformation 


(5.48) 


f^-h 

^auf 


xili = G 


(h) 

auj:,z 


+ --- + G 


(h) 

auj:,n 


^ (h) _ (h) 

oaup ,2 oauj:,i—1 


with 1 < j < n and gauf.o := 0- We regard as a transformation on all the 

coordinates, extending it as the identity on the remaining ones. We denote the new 
coordinates as 




^auf 

(h) 

Yauj 

x^G 

-^aur 


... 99(^1 ... 99(^1 1 

V^aufjlJ J ^our,n—ft —1’ ‘^aur,!’ ’ ^auf,n—h—l/ 

n>4, 2 < h < n — 2 
0 otherwise 

(„('*) ... „('*) ... „('*) ) 

v-^auj:,n—/i’ ’ >^auj:,n—1’ ^auj:,n—/i’ ’ ^auj:,n —1/ 

p,(ft) .(ft) .(ft) \ 

Wauf.li ■■■ I '-'aujr.m ^'-aur,!’ ^ 

. . . /(G \ 

l&auj:,!! j Bauf.n, «au!:,li J •^our.nj 


the new Hamiltonian as 


(5.49) 


r /i.(^) W W (h) \ n j^Yi — h/Ah) ih) (h) (h) \ 

/5ec,int,auy,/i—1 (taufi Yauj:, Xduy) :— /sec,int,/i —1 ^ Yaxi^ (^auy? ^auj:? Yau];? ^auy)- 

Now we define the domain where we want to consider /sec,int,auy,^-i- Firstly, we let 

—hAh.) (flG (h) (ft) \ 

^int.auy •— J auy? 5 J* 5 -^* 7 ^ ^int j 

where is defined in (5.23). Then ©(“nt.aur given by 

Sfnuauf = X X X (G* X X X 

with 


(G*)ci.e,Chr ■— (Gl)c 9 ejX---X {Gn-h)cn-hSn-h X {^n-h+l)ln-h+irn-h+l X ‘ ‘ ‘ X (G„_]^) J, 


-ITn 



5.2. SECULAR NORMALIZATIONS 


57 


Next, for 1 < h' < h and any fixed K > 0 and r > 2, we define 


(5.50) 


.n-h' (Ah) N 

\y auf,n—h' J ' 

+ p(l) p(l) 

2 ,n — 


^-IrAh) 


^, ('>') n 2 , ^ ('i') ^2 

+<%u,.n-fe'> pC/.') 




^scc (yauy,n —1 

IJi-h'(Ah) . 
^^scc \y avLi,n—h) 


n > 3, h' = 1, 2 < h < n 


n > 3, 2<h'<h<n — 

{h', h) ^ {n — l^n — 1) 
h' = h = n — 1. 


We then choose the following sub-domain of 


int, auy 


D 


sec,auj: 


f (Ah) (h) (h) (h) . 25 

\ V^Tiorm,auy? ^norm.auy? j norm,auj:?-^norm,auj:^ t 

7 


h 

int.auy 


n—h' 


k|> — 


(5.51) 


V k G \ {0}, \k\i<K, \/2<h' < h}. 


Here j is chosen to be 3, 2 or 1 accordingly to the three cases above. The set 
®int auf i® non-empty, if 7 is chosen suitably small. Indeed, if we put 


Ah) 

y auT,n—h 




(„A) 'la I (Jh) X 
\yaut,n—h/ ' \^aux.n—h/ 


aux,n—hJ p(^) 


(^(h) Ah) Ah) \ 

’ '^auj:,n—/i’ n—h^ ^ auj:,n—/i+l J 


then standard quantitative arguments show that, for any fixed value 


(\A) Ah) i e n 

w''auj:,n—/i’ ^ auj:,n—/i+l/ 


j.{h) ^{h) ^int.auj:; 

aup,n — h’A aup,n — h+l 


the measure of the set N„_,, C x Gn-h of 

where the inequality in (5.51) does not hold may be bounded as 


7 / 2 

measNji_/i < — meas [B- 


Gn-h), 


(where c depends only on the semi-axes ratio and the masses), hence (5.12) follows. 
This is because (yitr n-h) i® ^ diffeomorphism (Compare Appendix C). 
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Now we inspect the form of /scc,tnt,aur,/i-i in (5.49). Introducing the following 
symbols 


(juy 


1 

’ ^auj:,n — h—l’ ^auj:, 2 ’ ’ ^auj:,n—h—l/ 


(h) 

1 

n>4. 2 < h < n — 2. 

l<i<n — h — 1 

'aux,i 


otherwise 



.(h) ,_ ( 

auj:,z • \ 

'p w 

>'^auj:,25 

p(/i) a(^) 

' ' ' 5 '^aux,nf -^^auy,!’ 

... ) 

.(0 ._ / 

auj :,2 • \ 

.oauj:,2’ 

(h) 

■ ■ ■ ) gauj,nj 

Ah) 

•''ouf,!, 

pW 

7 ^aux,n) 

.(h) . , 

(P^"^ ■ 

' V/^auy,z 

f 

{P 

, /Ah) n2 

auj:,n—1/ ' Vyauj:,n—1/ 

auj:,z ’ \ 


2 

1 1 

2 


p(^) p(^) 




a(^) ... ) 

^'■auy,25 5 ^'-auy,n; 

() 7 r> 4 &l</i—l<n — 

v&auf.i+l Snur.i’ > &auf.n-/t &aur,n-/t-l/ u ^ ^ Cii 1 ^ u 


otherwise 

f(^) r^i^) 


X 




(h) 


(Cm, •••,Cu;.n-J n>4&l<h-l<n-3 


otherwise 

-.W 


~‘au^,n—h 

~(h) 

y 


(/i) _ (h) ■ (h) 

notm,_7 = Pn ormj ^^novxnj 

A2 , ( iri) .2 fjh) ^2 _L fjf^) 

y^aux,n—h-\-lJ ' \^aux,n—h-\-lJ 


V"auf,n—?!, ^auf,n—hji 

n('") '\2 , ( (h) 

aux,n—h-\-l/ ' \^aux,n—h-\-l 


(jyW )2 , / (C \2 

\^aux,n—l/ ' \^aux,n—l/ 


auj:,z 

— V 

2 

7 7 

2 


p (^) 

'^auy,2’ 

p(/i) 

* * • 5 '^auy,n7 ^'-auy,2’ 

^{h) \ 

7 ^^aux,n) 


Ah) 

Youp 

- 

Jauy.l’ 

-^auy •— -^auy,!’ ^^aup 

_ 

^auy,!’ 



by means of (5.29), we have 
(5.52) 

f (Ah) (h) (h) (Ah)^,r (Ah) (h) -{h) .{h)s 

/scc,int,auj:,/i—iCduy? Zquj:? Yauj:? ^aujrj —^sec,/il,y ) ' Jnorm,int,auj:,/i —iCduy^ y^up, X^ufj 

n—1 

= ^«c(C) 


i—n—h 
n—h— 1 

+ ^ /n 

2=1 


norm,int,auj:,/4 —1 (tauf.* ’ ^aux-i'n—h'! y auj:,z’ “^auy,^ 


(h) (h) 


~ih) ^(h) , 

— j aur,z’ ^aur,2/ 


in — h 


where we have let 

(5.53) 

/norm,ini,our,/i-l /norm,ini,/i-l O'/ouj , /norm,int,aur,li-l /n 

On the domain S 5 ec,au{: specified in (5.51), we aim to construct and real-analytic 
and canonical transformation 

(5.54) 

An-h . /Ah) (h) (h) (h) ] ^ J)h _i- It^^^ Ah) (h) Ah) \ ^ j)h 

rnorm.ouj ■ Vinorm,aujMnorm,aujj jnorm,aujMnorm,aup7 t ^ O^nuj) ^oujy nupMour^ ^ '‘-^inl 

such that the transformed Hamiltonian 

x—h 


fscc^aux.h •— /sec,int,auj:,/i—1 ' 


‘^norm,auy 
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has the form 

r _ 1 /-'(/l) \ I r ^ ^ 

/scc,auj:,/i — ^5zc,hyYnorm,au^) ' /norm,auj:,/i'^tnorm,aup5 ynorm,auy5 ^norm,auj:j 


n—1 


n—h— 1 


= V h* ') + V r ) 

/ j ^^scc w norm,auj:, 2 / ~ / ^ ‘/notm,auy,/iV^notm,auj:,i’‘ynotm,auj:, 2 ’^notm,auj:, 2 / 
i—n—h i—1 

r (^) (^) (^) ^ 

"T /cyp,see,auj:,/ii^notm,auj:5 Z^otTTx,auj:5 ynotm,auj:5 Xnorm,auj:j 


where 


' _ g-('') . . . _ rr(^) 

\onorm,auy,i+l o 


norm,aup ,2 


■ norm,aup ,2 




- ( 


notm,auy,i’ ’’’ ’ §norm,auj:,n—/i—1 ^novm,au^,n—h—2j 

qquadii n>4 & l<h—l<n — 3 


otherwise 


(Ah) '>2 , ( W .2 

v^notm.aur.n—/i/ wnotm.aur.n —/i/ 


, (• ('«) ',2 

\r'norm,auy,n—1/ ' \ynotm,auj:,n—1/ 


2 ’ ’ 2 

(/t) a(/i) »(/l) 


rW rA) \ih) xw N 

'^notm,nuj:,i> ’ '^norm,aur,n, ■■■ , Jrnorm.auy.n j 

and /cjp,scc,our,/i satisfies the bound for /cjp,scc,/i in (5.20). This will conclude the 
proof, up to apply the inverse transformation of (5.48), with 
replaced by and to take 

We shall obtain the transformation 4>nAm out i^ (5.54) via an application of Propo¬ 
sition D.l. Before doing it, we just remark that, since, in our particular case, 
./norm,int,auf./i-i depends On zitt yiul, xiul only via yitt and is even 

in (tatt n-h)’ proof of Proposition D.l can be easily handled to show that 
^nmm.auf oan be chosen of the form 


in—h 

^norm,aup • 


@('*) . = 
au^,J 

- - TT = 

oiiyj 


p(/t) (Ah) (h) ~ih) -(h) \ 

^ notm,aupj pnorm.ourj ynorm.our, Xnorm.aup) 

p(h) (AG (h) ~{h) -(h) \ 

'^notm,aur,iG>v>™.0“J' ^notm,auf,ri-h’ ynfrm.auf, Xnotm.ourJ 

j = 1, • • • , n-h-l 


AG 


ryyn) V-/ V -/ V -/ 

^notm,aupitnotm,auj:5 ^norm,auj:,n—/i’ ynotm,auy5 ^norm.auyj 

(X'G -(h) \_ -^(h) / (h) (h) -(h) -(h) X 

V^auy’ ^aup/ ~ ^norm,auyV^norm,auy7 Jnorm,auy5 ^nocm,auy; 


.{h) 

^norm,j 


AG 

^ auj:,n—/i+l 


/,(M (^{h) 

j = n — h + •••,n—1 


Xh) ft''*) 7<'*> v('‘> xf**) 1 

7'- ' p^notm,auj:,j''notm,au}:’notm,auj:,7i —/i’‘ynorm,auj:’ notm,auj:/ 

norm,auj:,j?^ 


X 


{h) 

auj:,n—/i+l 


AG 

A norm,auj:,n —/i+l 


JG 

Na 0 rm, auy ,n—1 


“rV^notm,auj:l^notm,aup5 ^norm.aup,n—h-’ ynotm,auj:5 ^norm.aupj 
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wViprpF*-^^ anri arp nrlrl- 7//^^ anr\ 

wiicic r notm,auj :5 'Jnorm,auj: ci-iiu ^riotTn,aup uuu, -A.riotm,auj:5 r notm,auj:,j ‘f^notT 

are even under the change 


(Ah) (h) 

V^norm,auj:5 ^notm,au^,n—h 


y (tnfirm.aiiri Z 


A) (h) . 

notm,auy5 ■^norm,auj:,n—/i/* 


Then (5.32)-(5.34) follow. 

Now we proceed with proving the existence of (pnatm auf We can choose, in (D.1),(D.3) 
and (D.4), 

(5.55) 

Pj = 2 (/i + l), ii = h, 
n — 1 

h(p,g,/)= h;„(yf^), 


rrii = 3i, i=l, — h—1 = N 

n—h— 1 


i—n—h 


2=1 


fHv n o 09 I •= p-'*-* (fA) ^{n-) ') 

J \ iiF'iHit) ' norTn,int.auj:,/i —1 ^ ^auj:,n—/i —2 ’auj:,n—/i —2 ’ ■^auj:,n—/i— 2 / 

3 := 3* := {(fc', k\ k'") &lAx IA+^ X = • • • = k'^_^ = 0 

= ... = fc" = 0 , C-. = ••• = C = 0 , k'i + ••• + = o} 


Ah) 


Ah) 


£ := 


:',fc",fc'")e3: fc;_, = fc"_, = o} 


where we have re-named 


(P>9) 

I 

Ui 


(Ah) Ah)^(Ah) ___ (h) (h) ___ (ft) X 

V/^auf? yauj:/ v^auj:,n—/i’ ’.^aup,n —1’ ^auj:,n—/i’ —1’/ 

/p(^) p(^) A A ^ 

y^au^,n-h^ *'* ’ , ii-auj:,n; 

r (^i) ,.. Ah) Jh) ^ ^ _ Jh) N 

\oauj:,n—/ 2 ’ ’ oduy,?!? 5 

{I, VAC), v-=v\ ^■■=e 


with 


77 := 


f := 


(G)auy,n —/i—2j ‘ ’ 5 G)auy,n —1; Gauy,n—/i—zj ' ‘ ’ j Gauj:,n—^ —Ij 
-^auy,n—/i—2; ‘ ‘ ’ 5 /i—l) 

(^auj:,n—/i —25 * * * ; ^auj:,n—17 §auj:,n—/i—27 * ‘ * 5 §auj:,n—/i — 1 ; 
^au^,n—h—it ‘**5 '^auj:,n—/i—l) ■ 


In order to verify that Proposition D.l can be applied, we have to check conditions 
(D.2) and (D.7). Due to the choices of 3? ^ and to the fact that only the function 
h"-'* in the summand for h^ec in (5.52) depends on (p^.n-ft^ ^^.n-ft^ G|,^tn-ft)’ 
it is sufficient to check that condition (D.2) holds with 

^ = <7^, (fc',fc) G Z2\{0}, K = K. 

But due to the choice of in (5.51), we have that (D.2) is verified, with 


a = 


«-ft)^ 


(®n-ft+l)^^"-^ ^ 


7 _ On-hA _ ra - 

f — ^ ~ ^hy dn-h- 


It remains to check the inequalities in (D.7). In view of the definition of /* following 
from the formulae (5.28), (5.53) and (5.55), of the definition of in (5.18), 
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the definition of the bound for in (5.20), and first inequality 

in (4.9), we see that the former of the inequalities in (D.7) is satisfied with 
(5.56) 


Ei = — max ■ 

Ch 




— } 


i = 1, 


n — h — 1. 


i^n-h-i+iy “1 ^n-h-i+1 

In order to check that also the second inequality in (D.7) is satisfied, we previously 
note that the number di in (D.6) can be taken to be 


J . r Qn—h'y Q 1 ... 

d^ = Chmin{^^^, Ori-h-i}, * = 1, 


, — /i — 1. 


Inserting then the above values for K, a, Ei and di into the left hand side of the 
second inequality in (D.7), we find that this can be bounded by 


1 

Ch 


■ max ■ 


(a+_^_,)2 (o„_^+i)3 j^f+i (a+_^_j2 (a„_^+i)3 


7^ 


{ay_hy (a„_,,_i+i)3’ 7 {a^-hy (a„_,,_,+i)3 

+ 3 

\ rr 


+ 
n- 
- . + 3 

—_ / n ' \ rr 




i—h 


7^ 


^n-hJ 


^2 n- 


'‘'n—/i—z+1 


On-h-i ■ 


K-h) 




Using (4.5), one easily finds that this quantity does not exceed 


(5.57) 


1 

— max 
Ch 


017^’ 7 


a' 


< 1 . 


where C/j depends only on the ratio and the masses and the inequality follows 

from (4.9). This conclude the proof of this case. □ 


5.2.4. Construction of T^ie arguments we have used in the previous 

section to construct ''') ^Iso fit for the case of therefore we 

shall not repeat them. We only limit to remark that, for this case. Equations (5.50), 
(5.55), (5.56) and (5.57) have to be replaced with 


d, 


,(i) 


\2 I /„(!) \2 

auj:,n —1-' ,n —1 

2 ’^auj:,n-l’^ciuj:,i 


"-i(y(^) D := 

c y-J auj:,n— 1 / 


h^ec(yilr.i) 


h"cc^(yilln-i) ^ > 3 


n = 2, 


f = /norm.inl.aur.o(tiuU’ yitf.U d, = Ci min { 


( 1 ) 


^-l7 


i = 1, • • • , n - 1, 6*0 := 01 


E,, = 


Cl 


max 


max 


{di) 


I 1 


an ’ 


{any 


{,.K{p 




i<-iy 


• max 




j^2t+2 


^2 


T 


\ 






i = 1 


n > 3, i = 2, 


n — 1 


□ 






















APPENDIX A 


Computing the domain of holomorphy 


A.l. On the analyticity of the solution of Kepler equation 

Here is a refinement of Proposition 4.1. 

Proposition A.l. Let e be as in (4.2). For any 0 < e < e there exists rj = rji^) 
such that, for any rj < rj < \ and any e G C with |e| < e, there exist two positive 
numbers ^ £ = £{r],e) such that the map 

(A.l) C G -f^(C)e) := C ~ esin^ 

is injective, its image verifies 

A:(T^-,e) VeeC: |e| < e. 

The inverse function 

£ gTj ^ C,{£,e) := K ^(^,e) G T^^(g) 

verifies 

(A.2) |1 — ecosC(^,e)| > 1 — ?7 

Therefore, C(£,e) is real-analytic for £ G Tj. 


The proof of Proposition A.l is elementary and goes along the same lines of [24]. 
Therefore, we shall present it skipping some detail. 

Lemma A. 2. Let e be as in Proposition f.l. For any 0 < e < e there exists a 
unique rj = rj(e) G (e, 1) such that 


V ?7 G [??, 1) : £r^{e) := log 


2 n 


- + \l 1 + Z2 


1 

e“ J 


-yrf + e^>0, f^(e) = 0 


77 = ry. 


Proof. By definition of e, and since the function p G [0,1] —?> 


creases with p, we have 


Consider now the function 


i+V^i+V 


1 + \/l + 


< 1 . 


p G (0,1] -)> gp{p) := 


pe 




77 + 

This function decreases with 77 for any p G (0,1]. Since 

gg\/i+®^ 


ge( 0 ) = e® > 1 , 5 >e(l) = 


1 + v'l + 


< 1 
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we find a unique rj = ri{e) € [0,1] such that 

9e{v) <1 V^<r?<l, g-eivie)) = 1- 

Since also 


5e(e) = 


e\/2 


= ^/2 


> 


we actually have 


1 + V2 1 + V2 

e <rj < 1. 


> 1 


□ 


Proof of Proposition A.l. We shall prove Proposition A.l with 


(A.3) 


C{r),e) := log 


£{'n,e) log 


+ ei + yjrf -e 


\/®i +' 


o2 


3 + i/1 + Z2 


1 

e‘J 


rf + 


where e = ei + ie2. Observe that i{r],e) > 0 by Lemma A. 2 . Moreover, since 

ei < |e| < e < ry < ry 

we have that C{ri,e) is well defined and positive^: 

C(?y,e) > log3 > 0. 
e 

We split Equation (A.l) into its real and imaginary part 

^1 (Cl, (2,61,62) := Cl - (61 sin Cl cosh C2 - 62 cos Ci sinh C2) = h 
Ar2(Ci,C2,6i,62) := C2 - (61 cos Cl sinhC2 + 62sinCi coshC2) = £2 

(with C = (1 + 1(2, £ = ^1+1^2)- The equation for the real part gives a unique 
solution 

Cl = Zi(ei,e2,C2,^i) 

provided 

(A. 4 ) |6i|<?y, |C2|<C(iy,6) 

since it reduces to an ordinary real Kepler equation 

Cl — Ai(ei,e2, (2) sin(Ci — C>i(ei,e2, (2) = ^1 if Ai(ei,02, (2) ^ 0 
(i = £1 otherwise 

with 

Ai(ei,e2,C2) := yjef cosh^ C2 + 6| sinh^ C2 

C'l (61,62, (2) : Ai cos (fi = ei cosh C2, Ei sin (f>i = 62 sinh C2 ■ 

and, under condition (A. 4 ), one has 
(A. 5 ) Ai < ry < 1. 


1 Actually, f(?7, e), as a function of (61,62), reaches its positive minimum 


Cmin = log [| + ^1+ > log(l + +2) 


for (61,62) = (0, e). 
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Observe that this solution Zi(ei,e2,C2,^i) verifies 

(A.6) Zi(ei,e2,-C2,^i) =-Zi(ei,e2,C2,^i) mod 27r. 

On the other hand, the function 

C2 K2(ei,62, C2,-^i) := Ar2(Zi(ei,e2,C2J-^1)1 C2,eij^2) 

is strictly increasing, therefore, it maps the interval [—C(^) e)) e)]; onto the 

interval [-L2 (17, e, £1), L2 (17, e, £1)], where 1^(77,e,£i) := K2(ei,e2,C(’7:e),-^i) 
(note that K2(ei,62, —C(i?,e),£1) = —K2(ei,e2, C(^ie),£1) because of (A.6)). We 
have thus proved that the map (A.l) maps bijectively the strip onto the set 

£ = £1 + i£2 G C : £1 S T, £2 S [—L2(?7,e,£1),L 2 ( 77 ,e,£1)]. 

But the curve 

£2 = L2(?7,e,£i) £iG[ 0,27r) 

is concave, its minimum points are cusps, where L2 attains the value 
L2.min(i?, e) = C(i?,e) - -el + el- 

The minimum of this quantity while |e| < e is just £{r],e) in (A.3). Inequality in 
(A.2) follows from 

|1 — e cos Cl > I Re (1 — ecosC)| > 1 — | Re (ecosC) | 

and (by (A. 5)) 

I Re (ecosC)l = |Ai(ei,e2,C2)cos(Ci - <('1(61,02,(2)1 < Ai < 77. 

□ 


A.2. Proof of Proposition 4.2 


Define 

Assume (4.3), with 
A:= (l-a^). 




6,:= \ll-el 


Q := 


^ (1 +<7)3(1 + <t 2)4’ 
Cl(<T)l5i 7 = 1, • • • ,77—1 


B := 


5n i = n 
s = ail — a) 


C,- := 


(1 -<72)(1+<7)3(1 +^2) 

C 2 (,a)\J^i +2g{ a)H] i = !,■■■ 
Si + ‘2g{cr)'^5l i = n 


,77—1 


where 

(A.7) 


Ci(<7) := v'1-<t2, C2(<7) := 
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and a, g are chosen as follows: g{a') is a suitable positive function, depending at 


most on the ratios such^ that 

A. ’ G.- 


(A.8) 


||c 0)||2 


g(cr') —0 as tr'—>■ 0, and [sinarg "^^^" | < g(o'')) J = 


provided 


max||arg(A,)|,|arg(xj)|,|arg(0j)|, |arg(i?j)|| < a' 


while a is so small that, ii £i, • • •, in are as in Proposition 4.1, with e replaced by 
ei, • • •, e„, then 


(T < min ■ 


and the following inequality is satisfied 
Cl (a) d, 


>1 V i = 1, • • • , n. 


Note that this inequality is satisfied for a suitably small, since, by definition, 
Sj > Sj, Ci((T')tl, C 2 (cr') i 1; as ct'—>■ 0. 

Definitions and assumptions in (4.3) imply, since a{l — a) < a, 

(1 - cr)G„ < |Xi| < G+(l + a) 

, , max I Im (Yj_i — Yi)| Oi 

tanarg(Yi_i — Xi)\ ^ -rr:—i-rr ^ - < cr < 1 

(A.9) |argx,|<|argx„_i|+ £ | sin"^ | < a < ^ 


i=j+i 


IXi -Xi+i| 


The previous inequalities imply that, firstly 

I I / 


< 


and, similarly. 


Xj-i (1 - cr)Gi; 


Xj 


< a 


therefore, the inequality for ij, ti is (4.4) follows. Secondly, the definitions of 0^, 
'df imply that conditions (A.14) are met and hence Lemma A.3 applies. By the 
thesis (A. 15), we have^, for j = 1, • • •, n — 1, 


^Since, for j = 1 ^ ||Cp^ |p depends only on Xj—i^ Xj-> ^rid 'dj as in ( 2 . 17 ) and all 

such coordinates, together also with Aj, have their anomalies bounded by cr', we can always find 
such a function p(cr'). 

^Beware that, if 2: = (21,22, zz) E C^, we denote 

\\z\\ .— Z-^ 22 “ 1 “ Z^. 

For a given 2 E C, the symbol \z\ denotes the usual modulus of 2 E C: 

1^1 := y(Re2)2 + (ImDX 
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|lic«f| 

< 

^2 + (1 + <^)(1 + ) Xi-i Xi 


< 




C, C,B2 

(A.IO) 

< 


For j = n, 






We suddenly have the left bound in (4.4): 


1 - 


"i.p 


< |1-' 


"i,P I 


|||C; 




A? 


< (5,- = 1 - 


for 7 = 1, • • •, n. Now we check the right bound. To this end, previously check the 
following inequality 


(A.ll) 


|lx,-i|-lx,l| > 


Because of the second inequality in (A.9), 

I arg [iXj-i - Xj){Xm-i - Am)] I < 2tan-i ct. 


Then we have 

1 _ ^2 

R-e [{Xj-l - Xj)(Xm-l -Xm)] > ^2 lAi-1 - XjllXm-1 “ Xm I • 

Taking the sum for m = •••, n, gives 


Re(Xj-i-Xj)Xi > 
> 


l-a2 

1 + Cr2 


IXj-1 


n 

-Xj\ lAm-l 


l-a^ 

1 + Cr2 


GJ\X,\ 


Xm I — 


1-S2 

1 + S2 


IXj-1 


XillXjl 


So, Lemma A.4 with 


^ = Xi-i, B = Xj, A = Gj-, 


a = 


l-a^ 
1 + 0-2 


gives (A.ll). Then the thesis (A.16) of Lemma A.3 and the definition of dj provide, 
for j = 1, • ••, n - 1, 


(A.12) |||C 




> 


A2cn 


A2(G-)2-(G+)2|i1,-^P> 


> {5] + V2g{a)5.){K+f 


where g{a) is as in (A.8). Again, this inequality is implied by the definition of 
19^ in (4.3) and the ones of A and C 2 in (A.7). By (A.8), (A.10) and (A.12), for 
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j = 1, • • • , n, we have 


ej.P = 


< 


(A.13) 


\ 


(1_ Re^^^)2 + (im 








A? 


< \/(1+ 2i5^3(cr)2 < 1 - = e|. 


For j = n, 


|||c[r)f I = |x„-iP > {S. + V2g{a)5r.)(Ktf 


again implies (A.13) with j = n. 

The proof of the inequality on the right in (4.4) proceeds in a similar way. Indeed, 
starting with 


Kp|^= |i4*+i)f _24).4+i) + ii4)f 


> 




- 2 


U/p U/p 


Il 4 ^f 


and using (as it follows from Proposition A.l) 


4+^’IP = |a^+i(l - ei+i,pCOsCi+i)^| > (1-??i+i)^(a,+i)^ 


and analogous arguments as above to evaluate 
finds the ansatz. 


„w_b+i) 

U/p Jj-p 


and 


ii 4 ^ip 


, one easily 

□ 


Estimates. 

Lemma A. 3. Fix a number cr > 0. Assume that, for 1 < j < n — 1, 

(A.14) 

R-eXj(Xj-i-Xj) > 0, |0j| < crmin{|xj-i|, |xtl}, | Im (i?j - 7r)| < log(l + cr). 

Then 

(A.15) |||c4f I < + (1 + a)(l + a2)|x,-i||x,P, - 

0)ll2| ^ |lxj--il - IXjll" 


(A.16) 


1 + 0-2 


(l + cr)(l + cr2)|xj_i||XjPj - Tip 
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Proof. We use the formula (2.19). By Taylor’s, given a, b, z G C, with 
|z| < o-miutgfo,!] \a + t{b-a)\ 

yjb'^ — z'^ — \/ — z"^ 


< 


< 

We use this formula with b := Xj-ii o := Xj^ ^ with the observation that, 

for ReXj(xj_i — Xj) > 0; the function 

t e [0,1] \xj + t{xj-i - Xj)f = IXjP + “^tReXjiXj-i - Xj) + - Xjf 

reaches its minimum, given by min{|xj_ip, |xj P}) for t = 0 or t = 1. Developing 
also the function ra G C —>■ cosw around w = it, with g := w — n = gi + ig 2 and 
\P 2 \ < log(l + cr) 


+1(6 — a) 


f ^\J{a -\-t{b — a))^ — z'^dt 
Jo d,t 

ib-a) f 
Jo 

\b-a\ f- 


. + t{b — a)) — z'^ 


zdt 


\a + t{b — a)| 


zdt 


+ t{b — 


a - z 


\b-a\ 

\/l — cr^ 


cos la + 11 = 

< 


/■i (f 

I (l-t)^cos(7r + t(ia 
\\Q?e\^^\<\\g\\l + a) 



leP sup I cos(7r + £i')| 


and using again the second inequality in (A. 14), then inequality in (A. 15) follows. 
The inequality in (A. 16) is obtained via the second inequality in (A. 14) and 





\Xj-i-Xj\ 

llxj-ip - Ixjf I 

llxj-il - Ixill 

\/l + cr^ 


□ 


Lemma A. 4. If A, B G C and a, A G K+ verify |A— il| > A and He B{A—B) > 
a|i?|A, where 0 < a < 1, then ||A| — |il|| > oA. 

Proof. Let D ■= A — B. Then | |A| — |i3| | = | |i? + D| — |i3| | < aA implies 
+ |D|2 + 2 Re'BD =\B + D\^ < {\B\ + aAf = \Bf + a'^iAf + 2a\B\A. 
This contradicts assumptions \D\ > A > aA and KeBD > a|i3|A. □ 




























APPENDIX B 


Proof of Lemma 3.2 


In this section, we prove the formulae (3.8) and (3.9) given in Lemma 3.2. 

We recall the following result 

Proposition B.l ([ 28 ]). Let X = Xi x • • • x c K® x • • • x K® and let 

(4,Xfc)eTixXfe^(4''^(4,Xfc),4''\4,Xfc))eR3xK3 k = i, ■■■, n 
be mappings such that, for 1 < i < j < n 

(A) the map 

(j)ij : {ti , 4, X„ Xj) {yf , yf> ,x^^\ xf') 

is symplectomorphism o/T^ x Xi x Xj into R.^^. 

(B) The map i£j,Xj) {yf\£j,Xj),x^^\ij,Xj)) verifies 

\\yf\e„x,)r _ ^ _n^. 

||xf (4,X,)|| 2A2 ’ 

where Aj is the variable conjugated to 4 ™ this symplectomorphism. 
Then the function 

P«(4,X) := -^Ide, 

3(4^^(£.,X0 -3:^(4-,X,))^ - |4g^(£.,X0||^||a:y(4-,X4|p 

2|l4'^(4,X4r 

is given by 

p . p« ^ aiT,-m2 3(,g).cy)^-||4-)p||cyf 1 r 

4 ||c«|14 2^yT||x«p- 

with C«(X) := xfie„X) X 4^'\4,x). 

Even though the (i,j) projections of the P-map do not verify assumption (A), one 
has 

Corollary B.2. The formula (B.l) applies also to the P-map, or, more in 
general, to any Kepler map K related to the map Pet in Definition 2.14 via 

XDe^ =1?(X). 

Proof. De£ verifies (A) and (B). □ 
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In particular, we have an expression for the second-order term of the doubly 
aged Newtonian potential 


— 2 
/k 


Mid£i 


miTTij 
(27r)2 jj2 


2||x«(£„XK)r 

Corollary B.3. For any K as in Corollary B.2, 
^2 


^(2) 




(B.2) 


= miiTij- 


+ i(5-4 


A3 


lie 


K 


,5 3||Cy|p, (j) 2 
^2 2^^^" ^ " 




^2 ; V- K 

Proof. Lemma B.2 implies that 

(2) Mim? ^ It (3(2^k^ 


(-i(i)\2 I 3||Cy||3 (j) 


c“f 


/] 


miruj 

1 

X — 
27r 




M, 




Mj 


/T||x«p- 


By (2.1) 


Ad cU) 


= 4.kp4+ 

= s^^,KP4 ■ ■ c\p 


Therefore, squaring, .^j-averaging and using 


we obtain 


[ia,^Kfd£i 
Ztt Jj 

= 


F [(hi^K)‘^d£i 

ZTT Jy 


Pl|c«p 


2 Af 

ZTT Jj 

= 

0 


of 

2 


-f 

«? 



11^^ 

A? ' 


'-)(p 4 


U )\2 


C^i>) 


Using finally 


Kk pdU = - 2 


5 3 IIC^^p 

A? 


h 


d£-i 


'T X 


A) 112 


1 A, 
“'lic^i 


aver- 


we obtain (B.2). 


□ 
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Now we may proceed with proving the formulae in (3.8) and (3.9). 

Proof of of (3.8). We apply Corollary B.3 with K = P, i = n — 1, j = n. 
Using ||Cp"'^|| = Xrt-i (see (2.17)), Cp"^ = Sp”^ and Eq. (2.3), Proposition 2.6, and 
Remark 2.7, we have 


(n-l) q(u.) 


Pi 


j(") 


0n-l 

1 


(n-l) 






g(n-l) 


-S^”0 X 
(n-l) _ g(n) 


)).s[ 


X P 


lie 


(n-l). 


Ui - Ql-i)ixl-2 - ^l-i) sini9„-i. 


Proof of (3.9). By Corollary B.3 with K = P, j = i + 1, we find, for /j 
an expression as in (B.2), replacing (n — l,n) with (i,i + 1). 


2 , 2+1 


□ 

( 2 ) 


(B.3) 


pW , (^d+i) ^ pW , ^gd+i) _ gd+2)^ = 0, - p, 


Q 


{i) p(2+l) 


= • (s 


p 

(*+i) 


-s 


p / 
(*+ 2 )', 


(i) q(i+2) 

P ■ ^P 


||C«| 


?-0f)(xLi-0f)sini?. 


c(*) w pco pco ,, qC'-i-j-; 

Op A Jr p • Op Jr p A Op • Op ). 


j{i) q(i+2) 


jCd 

P 


;(i+l) 

^P 


Now, when (0i+i,i?i+i) = (0,7r), ||Cp''"^^|| reduces to 




Xi+l; 


(provided arg(xi - Xi+i) G (“f, f] mod 27r) and 


so 


Mi+2) 


Xi+1 g(^+i) 
X* ^ 


and hence, the extra-terms in (B.3) reduce to 


p*) _ g^i+2) ^ Q 


X^+l 

X^ 


c(^) V c(^+2) 

Op A Jr p • Op 


Xi+i 
Xi 

X = o. 


xLi-0f\/xf-0fsini?, 


Then (3.9) readily follows. 


□ 















APPENDIX C 


Checking the non-degeneracy condition 


In this section we prove statement 4 of Proposition 5.3. 

Due to the form of h^ec in (5.I7)-(5.I8) and to the bound for h'^^. ^ in (5.20), it is 
sufficient to prove that the maps 


a 


w 


.A _o 


in (5.18), where 

(pW)2 + (^W)2 


cf ^ = 


■, x^) 

i = 1 

3 

II 

to 

, xi-2> xih) 

i = n — 1 

Sz n> 3 

xLi) 

i = 2 , • • • , n — 2 

& n > 4 


2=1 

& n > 3 




are diffeomorphisms, with non-vanishing Hessian matrices. We shall do this verifi¬ 
cations for just one of the cases above, and we choose the second case in the list, 
i = n — 1, for u > 3. The explicit expression of h"c7^ is given in (5.42)-(5.43). We 
neglect the coefficient A„_i (which does not depend on Ci-\) nnd we denote 


-un—1 T 7 , I o Pn—1 T 9n—1 , _ fPn—1 A Qn—1 \2 

-tlscc — -C-'n—1 “r iin— 1 ^ “r 


)^ -|- 0(p„_i, (7„_i)® 

z, z 

—rr f 

the function h"ec thus rescaled, and its gradient with respect to ( " ^ ^—: 
Xn- 2 , Xn-i)- A perturbative argument shows that, under the choices of Corollary 
4.3, the frequency-map with respect to (y„_ 2 ,Xn-i) associated to 


En-i — — 


A3 


5-3 


(Xn-2 - Xn-l)^ 


2Xn-l^ ^n-1 


is an injection of its domain and hence, by another perturbative argument, so is 

the gradient of h"^^^ with respect to the same coordinates, for any fixed value 
2 2 

of 2 ^^3 the other hand, since t„_i does not vanish under the same 

assumptions of Corollary 4.3, is an injection. The computation shows that 

the Jacobian of does not vanish. □ 
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APPENDIX D 


Some results from perturbation theory 


D.l. A multi-scale normal form theorem 


The purpose of this section is to present a normal form result which takes 
into account different scale lengths. It is a particularization of [31, Normal Form 
Lemma, p. 192] and uses the same techniques of that paper. 

Following [31], the notations are as follows. 

• If A C is open and connected, T := ]R/(27rZ) is the usual flat torus, 
r, s are positive numbers, we denote as := IJoigA ^ ^ ^ 

the complex r-neighborhood of A. will denote the complex set 

T + i[—s, sj. As usual, denotes the ball in <C^ with radius r centered 

at X, accordingly to a prefixed norm j ■ j of . 

. If / = f{u,p, q, if) is real-analytic for {u,p, q, p) S Wy^s,e = UyX x T^, 
and affords the Taylor-Fourier expansion 

j=i v2 iV2 

we denote as |[/|[„,s,e its “sup-(Taylor, Fourier) norm”: 

|[/|[„..,, := ^ sup 

with \k\ := \k\i, |(a,^)| := jaji + |/3|i. 

. If / is as in the previous item, K > 0 and £ = £i x £2 is a sub-lattice 
of X Z^, Tfcf and n^/ denote, respectively, the A-truncation and the 
£-projection of /: 


TKf:= 

\{a,l3)\<K, \k\<K j=l 


fceiii 




V2 


iv^ 


Proposition D.l (Multi-scale normal form). Let 

v, I, 1 < mi < • • • < niN = m 


be natural numbers; 

A c S c Cl, C[ C C 2 , C 2 C K™ 2 -rni^ ••• , Cat, C)v C 
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be open and connected sets; 

r, s, s, pi > P 2 ■ ■ ■ > PN, Pi > P2 ■ ■ ■ ^ p'n 
positive numbers. Put 

Vi := - - ,p'i), V := vn 


U — 

LJ V .— ^VN 


w — 

y* v,s,£ •— yyvN,s,£t 


C/« := A X Cl,, X ... X Cl,, X c;,, X ... X C',,, 

:= C/^f X Tr X B,, 
with 7 = 1, ..., N. 

Let a, K > 0 with 0 < s < 61og5/6 and Ks > 12; let also £, and 3i, • • •, 3 a be 
sub-lattices ofl} x and tei 3 := 3i U ... U 3 a. 

Let 

(D.l) H{u,ip,p,q) = h{p,q,I) + f{u,ip,p,q) 

be real-analytic for {u,p,p,q) G Wv,s,e, where u : = {I, p, = (Ji, ..., Iv, pi, 
■■ ■, Pm, fi, ■■ ■, fm). Suppose that 

2 I 2 

(i) h depends on {p, q) only via frequency map u) = (wi, ..., 

LOi, LUi+i, ■ ■ ■, uji+i,) defined via 


9p2+g2 h 1 < i < 


oji := 


dii_ih. £-\-l<i<£-\-v 

verifying 

(D.2) \io{p,q,I)-{k',k)\>a V (fc',fc) G 3\£, |(fc',A:)| < X 

and all {p,q,I) G x Ar,' 

(ii) / is a sum 

N 

(D.3) f = '^Mu„p,p,q) 

where fi is real-analytie on Wvl]s,e cind has the form 

(D.4) 

Mu„p,p,q) = Yi n 


(a-a+,fc)e3i 


1=1 


fc=l 


V2 


V2i 


with 


(D.5) Ui := {I, p\ C) ■= ih, ■■■ , 4, Vi, ■■■ , Vmi, 6, ••• . ?mi); 

(iii) the following “smallness” conditions hold. If 
(D.6) Ci := e(l + Ce + mje)/2, dj := min{rs, e^, pip'} 

with e denoting Neper number, then 


N 


(D.7) 


ll/.llvu« < r 






adt 


»■+ 
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Then, one can find a real-analytic and symplectic transformation 

‘i’ • W^D/6"',s/6^,e/6'V Wv,<T,e 

which conjugates H to 


N 


N 


H*{u,(p,p,q) := iJ o$ = h{I,p,q) + '^gi{ui,ip,p,q) + '^ f*{u,(p,p,q), 

i=l i=l 

where gi, fi verify 

9i = n3.n£Tff5i 
Il5i - n3,n£'7A/i||„,/6N,o-/6'v,e/6N < 

, \ ^ ||/j|Uj,s,e II J |, 

^ adj 

j=i j 

‘^‘^C.k\\fk\\vk,s,eKS 


^8'^ adk 


I .s.e 


k=l 


\\f*\\vi/6r^,s/6^.e/6^ < (g)^ \\fi\\vi,s,e 

Finally, <i> is close to the identity in the following sense. Given F, real-analytic on 


W 


(*) 


Vi/6’^,s/6^’ 


||i^O$-i^||„/giV_s/6N-,£/6iv < ^- k\\fk\\vk,s,e - \\F\\.„./qN ,s/6^ 

k—1 ^ ^^k,i 

with dk,i ■= Taax{dk,di}. 

The proof of Proposition D.l is based on the following 

Lemma D.2. Let N gN, v, i, rm. A, B, C^, C[, r, s, pi, p[, Uvf , Wvl]s,e, Ci, 
di, with i = 1, ■ ■ ■ , N+1, be as in Proposition D.l; v := {r, pi, ■ ■ ■ , PA+i) Po ' ‘ ‘ ) Pa-i-i)’ 
U — ur — Let 

(D.8) H (p, q, I, ip, p, = h(p, q, I) + g{p, q, I,p,p,f) + f{p, q, I, p, p, 

be real-analytic for (u,p,p,q) G Wy^s,e- Suppose assumption (i) of Proposition D.l 
and, moreover, the following ones 

(ii) g is a sum 

N 

(D-9) 9 = '^gz{ui,‘p,p,q) 

where gi is real-analytic on Wvl]s,s o-nd m is as in (D.5); 

(iii) 91 , ■■■, 9n and f satisfy 

9i = HtPj, / = n3/ 

and 


N 


2'^c^Ks , 


(D-10) E <!’ 


2 = 1 


< 


ad 


iV+l 


^ i^cpk+iKs 
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Then, one can find a real-analytic and symplectic transformation 

$ : (m', (p',p', q') e W„iQ^s/e,EtG («, q) € Wv,a,e 

such that 

■= Tlo^ = h-\-g-\-g,ic-\- f^, 


where g* = Ti^risfTKg* *s 3 H Z-resonant and the following bounds hold 

,2^Cjq_^l ll/llw.Sje , 2’^Ci ll^illi 


II 5 * ^K^jr\2f\\v/6,a-/6,e/6 — (" 


ad 


A+l 


E 

1=1 


adi 


< 


|t»,s,£: 


ll/*IU/6,cr/6.E/6 < e 




Finally, $ is close to the identity in the following sense: for any F which is real- 
analytic on Wv'^}s^e, 


(D.ll) ||F O $ — $||„/6,s/6,e/6 ^ 


‘^'^Cf!+f\\f\\y^s,eKs 

adi 


ll^lUi.s.e < „ ||^|U,s,e- 


The following Lemma is a trivial extension^ of [31, Iterative Lemma]. Its proof is 
omitted. 

Lemma D.3. Let s = (si,--- ,s,j), r = (ri,--- ,r,y), e = (ei,--- ,ei), p = 
{pi,--- ,PTn), p' = (p'l,--- ,Pm); '■= (.r,p,p'), V := {r, p, p') < v/2, s < s/2, 
e < e/2, 

fe = l,--- ,m 

Let 

m 

H{u,ip,p,q) = h{I,p,q)+g{u,ip,p,q)+f{u,ip,p,q) g{u,ip,p,q) = ^gi{u,ip,p,q) 
he real-analytic on Assume that inequality (D.2) and 

(D.12) ||/|U,s,e < - 

c 

are satisfied. Then one can find a real-analytic and symplectic transformation 

$ ■ ^^v — 2v,s—2s^€—2i ^ ^^v,s,e 


1 


In order to obtain the extension it is sufficient to replace (p of [31, Appendix A] with 


<p= E 

(cx-0,k)G^\£i 
\{c.,0)\<K, \k\<K 


fk,a,p{u) TT , Pj - ^1j -sa, f Pj + ^1] -.g, 

i{a- 0 ,k)-u j-r ^ ^ > 
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defined by the time-one floufi X^f := / o $ of a suitable verifying 

IIJ.II ^ ll/l|j;,s,e 

a 

such that 

H+ :=Ho^ = h-\-g-\- Il^ns/ + /+ 
and, moreover, the following bounds hold 


II/+IU— 2fi,s—2s,e—2e ^ (l 


aJ 


LJ,s,e) 


-O^ 


+e 


-Ks 


v,s,e ( ~ ) 11/11?;, s,e + II {'5^; ff } II?; —?),s —s,e—e 


Finally, for any real-analytic function F on Wv,s,e, 

l|{0; -^111?; —?;,s —s,e—e 


|i^OlI> F^y—2v,s — 2s,€—2i — 


1 - 


0(5 


Proof of Lemma D.2. Following [31], the proof is obtained via iterate 
applications of Lemma D.3. 

To avoid too many indices, we shall prove this lemma taking, in (D.9), 5V = 1; the 
extension to 5V > 1 being straightforward. Namely, we take 

Rl — ' ' ' — Rmi = P, Pi — ' ' ' = Pmi ~ P 
(D.13) Pmi-l-l = • • • = Pm = P, Pm^ + l = • • • = Rm = P 

where 1 < mi < m. Letting 

v:={r,p,p'), v:={r,p,p'), E:=\\f\\y^s,s, G := \\g\\y^s,e, c=ci, 

J := min{rs, e^,pp'}, d:={rs, e^,pp'}, 

we rewrite the assumptions in (D.IO) as 


(D.14) 


2^cGKs 


< 1 , 


2'^cEKs 


< 1 . 


ad ’ ad 

The inequality on the right clearly implies (D.12). So, we apply Lemma D.3 to the 
Hamiltonian (D.8), taking n = ■ ■ ■ = r,^ = r, si = ■ ■ ■ = Si, = s, ei = ■ ■ ■ = ee = e, 
Pk, p'j^. as in (D.13) and 

V = Vo := v/6, s = So := s/6, e = sq ■= e/6 

V = Vo ■= il/6, s := So := s/6, e := Cq := e/6 

6 := {fs, e^,pp} = 

Letting 

vi := V — 2vo = 3/4z;, si := s — 2s = 2/3s, ei := e — 2i = 2/3e 

by Lemma D.3, we find a canonical transformation <i>o = X,pg which is real-analytic 
on Wy„^si,ei and conjugates F[ to Hi = h + g + gi + fi, where gi = H^nsPif / and 


^The time-one flow generated by (p is defined as the differential operator 




k=0 


k\ 


where L9/ := / and L^/ := {0,L^ ^/}) with k = 1, 2, 
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II/: 




< ( 1 - 


36cE 

ad 


,„-Ks/6 , 

V ad ^ 


E 


where 


Here, we have used 


36c£;,_i36cG^ 

(1 — -:— ) — E 

ad ad 


6 := inin{fs, i^,pp'} = 


36' 


ll {?^5 9 } i ,^ ,i^^^\\v-v.s-S,e-i — 3} j ||c_5_s-s,£-e 


(D.15) 


cG cG 


since g depends on 77 , ^ only via 77 ^ = {pi,--- , 77 ^ 1 ), = (fi,--- ,Cmi)- It is 

sufficient to consider the case 

-Ks/6 , 75 ,^ 18cE 
- ad 

since otherwise the Lemma is proved. In such case, using (D.14) we can write 
^ ^ 32. 9 2'^cEKs , 9 2'^cEKs , 9 2'^cGKs^ E 

-1^1 — \\jl\\vi,si,ei S 


(D.16) 

Let 

Note that 
(D.17) 


ad 64 ad 32 ad Ks 

E f 2’^cEKs 2'^cGKs '1 E 

< max <-, -=— > < — 

Ks I ad ’ ad / 4 


L ■= 


Ks 


■ 121 og 2 J 


L > 1, Ks > 8L, 


since we have assumed Ks > 12. We want to prove that Lemma D.3 can be applied 
L times with parameters 


(D.18) 

For L = 1, this follows from (D.16): 


4L’ 4L’ 4L’ 


d,. = 


16L2^ 


i = 1 , 


.L. 


E. := IIMI,,,..... < 


< 2 -“— 

C 


which is implied by the inequality in (D.16) and assumption (D.IO). We then 
assume L > 2. Suppose, by induction, that, for a certain 1 < i < L — 1, and any 
1 < J < 7 , we have conjugated H to 


Hj —h + g + gj + fj 


where gj = XilJo H^nsTif/fe 


(D.19) 


^7 := \\f 3 \\vj,s„ej < 2 
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where vq, sq, £o are as above, vq := v, sq := s, eq := e and Vj = Vj-i — 2vj-i. Then 
by Lemma D.3, on the domain we fined a real-analytic transforma¬ 

tion ihi = X^., which conjugates Hi to 


Hi+i — h -I- 5 -I- gi+i + 

where g^+i = gi + Hsmaft = Sfc=o n£n3^A/fc- We prove that (D.19) is satisfied 
for j = i + 1. Using^ the assumption on the right in (D.14), (D.16), the inequality 
for Ks in (D.17) and the definition of Si in (D.18), we have 


\\{9i^4’i}\\vi-Vi,Si-Si,ei-ei < (-E"! + -^)]-E'i < 


C E C E 

-a(5i Ks aJ,; L. 


E, < 


El 

32' 


Moreover, by a similar argument as in (D.15) and since g is actually real-analytic 
in the larger domain 




we have 




\Vi —Vi,Si — Si,Si —Si 


— 11/ A III <r 2EiE. ^ El 


where 


Si := min{fjSi,;0i/5-} = 


16L2’ 


i = I,--- ,L. 


Then we find"* 


cE- 

Ei+l — ll/i-l-1 ^ (1 ~ ) 


cEi 


0^1 1 a^i 




E, 


rE 

rE 

64.1 1 1 

- 63^64 8 32 64^ 


< 


E, 


,_6 


< E^i < 2 
4 c 


since i > 1. Then we let $ := $o o • • • o iL* :=Ho^ = h + g + gh+i + fh+i, 
9* ■= 9 l+ 1 j /* := Jl+i and we have, by telescopic inequalities and (D.16), 

L L L ^ 

\\g* -As,nfiTKf\\v/6,s/6,e/6 = E, W^SnSjTK fi\\ <E,Ei < EiE, TTET 


i=l 


i=l 


i=l 


- 3 ®‘- <Tr + 


^For the proof of inequality \\{gi, (l}i^\\y._{j.^Si—Si,ei—ei ^ t)’ f^ompare [ 31 , Proof 

of the Normal Form Lemma]. 

'^Since K > 8 L and L > 2 , one has (1 — ^ no- ^ I8n with the r.h.s bounded above by 

(4/7)16 (^^ decreases to as L —>■ +oo). 
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Now we prove (D.ll). Let F € Wg_s,e, := F, Fi := F o o ■ ■ ■ o i = 0, • • •, 
L. Then 


|i^O$ T’llu/e,s/6,e/6 — II^L ^lUi + i.SL+i.eL+i < O-Fi-1 ||t)i+i,s 


^ E 


2 = 0 

cEj c.E'i 

^_a£i^_|[p|l ^ Z^2^0 aSj 

^_c^Jl^llSLSi,e. - . 

i=0 '' aSi > aSi ) 


ll-^IU,s,e 


< f ^,ia.5||i.|i,„<?!£§L£||F|i,„ 


2=0 


ad 


where we have used < 1/24 that, for 0 < x < 1/24, log(l — x) ^ < |x and 


E 


i=0 


cEi 

aJi 


cEq cEi Qi^cEq cE\ 1 

Q^o ~ a^i 4®“^ 

2=1 2=1 

2^cEq 4 cL^i 2^cEoKs 

ad 3 a^i ad 


The proof for E G 144,s,e is similar. □ 

Proof of Proposition D.l. For simplicity of notations, we prove Proposi¬ 
tion D.l in the case u = .^ = 1; the generalization to any u, i being straightforward. 
Consider the Hamiltonian 


Ho{ui,(p,p,q) := h{I,p,q) + fi{ui,p,p,q), {ui,(p,p,q) G 
To this Hamiltonian let us apply Lemma D.2, with g = 0, so as to conjugate it to 
Hi := Ho 0^1 = h + gi + fil\ iui,p,p,q) G 
where 51 , correspond to (/,, /,, hence satisfy 



< 


,-i.s/ 6 ||^W|| 


Hi ,s,e 


llffi 


Then we have 


ll5l||Hi/6,s/6,e/6 — 

7 

g ll/l|Ui,s,e 

n£n3^-^/l 11 Hi/6, s/6,e/6 E 

2"ci||/i||2^,,,, 

adi 


E^\u,p,p,q) := iJ o $1 


N 

Hq o 4)1 + fj o $1 
/=2 


ii + 5i + /il 


N 

E/r 

i =2 


where := fj o d*!. Assume, inductively, that, for some 1 < z < N — 1 and any 
1 < J 4 * we have conjugated H to 


, (^,p, g) = i/ o $1 o • • • o 


3 

= h + ^ gt 

fc=i 


3 

E 

fc=i 


fU) 
3 fe* 


N 

E 

k=j+l 


fU) 

■Ik 


: W 


U) 


v/Q^ ,s/6) ,e/6) 


^^t;/6)-Ls/6)-),e/6)- 


where 










D.l. A MULTI-SCALE NORMAL FORM THEOREM 


85 


transforms 


into 


i-i 


Hj-i h Qk + / 


(i-i) 


k=l 


Hj-i — h + ^ gfe + . 

k=l 

The Hamiltonian 

i 

Hi [ui+i, If, p, g) := h + ^ gk (uk, f, p, q) + f-% [ut+i, f, p, q) 

k^l 

is real-analytic for {ui^i,(p,p,q) € s/ 6* e/6* satisfies the assumptions 

of Lemma D.2, with N = i. Then one can find $i+i : W 


W 


(*+i) 


Ui-|.i/6*,s/6*,c/6' 


such that Hi o $^+1 = h + ffk + where 


(i+1) 

Vi+i /6^+^,s/6^+^ ,e/6^+^ 


\\fL+l\\vi+i/6*+i-,s/6'+i-,s/6'+i^ < e \\fi+i\\vi+i/6',s/6\s/6' 


-As/6’+L| Ai) 




Hi + i,S,£ 


||5i-l-l||iii-i-i/6*+i,s/6*+i.e/6*+i ^ g II/i+l IUi-(-i/6Us/6Le/6* ^ 6^8^ ll/i-l-1 llni+ijS.e 

||5*-|-1 - n£n3'^-fS'/i+lLi-i-i/6*+i,s/6*+i,e/6*+i < IlSi+l “ n£n3'7R/i+\ Li+i/6*+i,s/6*+i .e/6*+i 

+ ||n£n3^-R'/f+i “ n£n3Tif/i+i||^._^^/6*+i.s/6*+i,e/6*+i 

< ||5,+i-n£n3rx/i;\|| t;i+i/6*+Ls/6*+Le/6*+i 

ll/i+l ~ /j+l|Ui-|-i/6*+i,s/6*+i,e/6‘+i 


9,0,- 2’’ci-|-i ||/i-|-l|| 


ViJ^l ,S,S 


ad, 


1+1 


< 


^3 WJ3WUj+,s 


3 = 1 


m 




adj 

yh 

adk 


||/i-l-l|Ui+i,s,e 


-ll/.+lll 


Vi^x,s,e 


with := fl^}o^ij^i for 1 < fc < i+1 and := o<^i+i for i + 2 < k < N. 

Then we find 


^d+i) H^i)oi^>,+, = {h + Y,9k + Y.fkI+ T. 


N 




i+1 


fc=l /c=l 


k—i+l 


N 


k—1 k—i-\-2 


i-\-l 2+1 N 

Ik* Jk 

k—1 k—1 k—i-\-2 


- ^+'^9k +'^/fc 
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(N) 


and hence, after N steps, 

hw H o ■ ■ ■ o = h+fffc + /j 

fc=i k=i 

satisfies the thesis of Proposition D.l. 

D.2. A slightly-perturbed integrable system 


N 


J+1 


□ 


The following result is well known in the literature of close-to be integrable 
systems, hence its proof is omitted. Note that it deals with an integrable system, 
close to another integrable one. 

Theorem D.4. One can find a number Cq such that, for any real-analytic, 
one-dimensional, system 


H(P,Q)=h(- 




) + f{P,Q) (P,Q) e® = i?2(0) cC^ 


and any 0 < e < e, such that 


(D.20) 


inf Idhl > a, 
B? 


sup I/I < e, 
B? 


Co 0£2 


< 1 , 


one can find a real-analytic transformation 

(/), : (P„ Q,) e ^ Q) g p2 

which conjugates H to a function H» = H o depending only on * . The as¬ 

sertion can be extended to the case that h, / are functions of other canonical coor- 
dinates (P', Q', y, x), depending on them only viaY = (y, ^ , ■ • ■ , ^ ), 

with y G Yp, {Pj,Qj) G P^/. In this case, letting (P*,(5*) —>■ (/*(P*,Q*;Y) t/ie 
transformation obtained for any fixed value of Y, there exists a canonical, real- 
analytic, transformation of the form 

: (P,Q) = (/)*(P*,(3*; Y*) y = y*, x = x*+(/j(Y*), Pj+iQ'j = 

which conjugates H to a function H* = H o depending only on * and Y,. 
In this case, the functions ipj, ipj verify 


\Tj\ < 


Co apj 


IV'il < 


1 


r ! 

Co ae' 








APPENDIX E 


More on the geometrical structure of the 
P-coordinates, compared to Deprit’s coordinates 


In this section we aim to point out differences and similarities between the 
P-coordinates and the coordinates denoted as ('h,P, A, 7 ,£) in [7, 27, 9]. 

We recall that the “planetary” coordinates ('h, P, A, ijj, 7 , £) may be derived (after a 
canonical transformation) from a more general set of canonical coordinates studied 
by A. Deprit. In their planetary form, the coordinates (4', P, A, ■)/), 7 ,1!) have been 
rediscovered^ by the author during her PhD, under the strong motivation of their 
application to the planetary problem [27, 9], 

Let us recall their definition^, in the spirit of Kepler maps (Definition 2.2). 

Let Cg\ be as in (2.8) of Section 2 and define the Dep-nodes 




(E.l) 


Hi := < 


c(0 c(*+l) 


^n —1 




i = 0 

j = 1, • • • , n — 1. 


Then let 


Eoep ■ — 


,•••,£„) C E^x-- 

■xE^) : 

0 < 

Ci < 1 , 

rii-i 7 ^ 0 Vt = 1, • 

On Ep)ep) define the map 






''Dep • ■ ■ ■ I 

G Eoep ~ 

^Xd, 

sp € ^Dep — 7'Qgp(EEep) 

where 









Xoep = («',r,A,' 


* X 1 

X Rl 

X T" X T” 

where 







= 

(^- 1 , 

tl'o, ^) G K+ X R+ 

X 

= 

(t/-_i,t/;o,§) G T X T X T"-2 

P = 

(Fi,-- 

• ,r„) gK!^ 


7 = 

(71G-- 

, 7 n)GT" 

A = 

(Ai, • • 

• ,A„) gR!^ 






with 


^' = (^'l,--- ,^'n_2) ^/> = (V’l,--- ,'0n-2) 


^The proof of their symplectic character found in [27] has been published in [7]. Another 
proof has been given in |36| . 

^For sake of uniformity, we use slightly different notations with respect to the ones in |7], 
actually closer to the ones of the paper |12]). 
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are defined as follows. The coordinates Kj are as in (2.11), while (4',r,'0,7) are 
defined as 

(E.2) 




f C := afe(3)(A:^^\no) 

i = 1 

II 

CM 

1 

1 

II 

1 



i 

1 

V Oj, 

2 < i < n 

jl 

10" 

7i := 

1 < i < n 


Then is a bijection [12, 27, 7, 36]. 

Definition E.l. We call Deprit’s map, or Dep map, the Kepler map 


Dep : Dep = (Xnep, G Doep = Xoep x T" ^ {y, x) G x 


associated to rQep- 


Comparing P and Dep 

a) Both the P and Dep-coordinates reduce the system to (3n—2) degrees of freedom. 
They share the following three coordinates (two actions and an angle) 


'S-i = Z = Qo, 


■0-1 = C = 1 ^ 0 , ^'o = G = Xo 


which are integrals of the system. As a consequence, the coordinates (^, C) and, 
respectively. 


g := 00, 0 := Ko 


do not appear into the Hamiltonian. Note that Dep and P share also the fixed node 
no = ui. 

b) The angle g for the set Dep describes the motion of the node ni in (E.l) and, 
by the cyclic character of g, this motion is negligible. Its counterpart in the set P 
is the node ni in (2.10), the negligible motion of which is governed by g. 

c) Compare the diagrams in (2.20) and (2.21) with the two ones associated to the 
Dep-map, respectively: 
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TIq TT-i . Tlfi—2 ^n—1 



fr 

tr 


tr tr 

fc(3) ^ 

s« - 


... ^ = 


i 

i 


i 


4 '’ 



E 



t 




-ni 

— 712 


^n—1 

and 





Fo 

Ft 

... 

^ F* 

^ ••• ^ F*=G* 




i 

; i 


Gt 


G* 


where 

F* = {ui- 

q(0 'i 

1 . ) 

G* = ( 

-n*, i = !,■■■ ,n. 

Note that, analogously to (2.20), n 
vectors in the tree (2.20). 

i in (E.l) is the skew-product of its two previous 

d) While Dep is not defined for the planar problem, P is, and, in that case, the 
coordinates (0 ,x,d, k) in (2.11) reduce to^ 

f Xo 



1 

f argP^F _1 . ^ Q 

0 .= < 

= 

= 

Ki = i 

1 ^ 

1 0 


1 ^ 

1 

[ prOpO+l) + TT 7 =1,-- -,77,-1 


x.= E iice^ii 

j=i+l 

while the (A, €) remain unchanged. 

e) The P-map is singular when some eccentricity vanishes or some of the following 
relations hold 

II p(0 II II 

^ Here by “planar case” we mean Cg ^ || ■ ■ ■ || || Note that, to be more precise, 

and K,o would not exist in that case (since ui =0). However, since they are both cyclic angles, we 
can fix them to an arbitrary value. The choice above corresponds to replace ui with 
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The former of such relations is negligible, while the other ones have no physical 
meaning. Therefore, the only physically relevant singularities of P are for zero- 
eccentric motions. 

The Dep-map is singular when some eccentricity ti vanishes or some of the following 
relations hold 

The configurations || have a relevant physical meaning, since the planar 

case corresponds to the intersection of all such configurations. A complete regular¬ 
ization of all the singularities of the Dep-map has been obtained in [27, 9], which 
allowed to overcome the problem of the rotational degeneracy (see [8] for informa¬ 
tion) of the planetary problem and to construct the Brkhoff normal form of it. It 
works at expenses of one extra-degree of freedom. 


f) The Euclidean lengths ||Cg^|| of the planets’ angular momenta are the actions T^ 
among Dep-coordinates: see (E.2). In terms of the P-coordinates they have more 
involved expressions in (2.17). As mentioned in the previous item, this makes more 
difficult regularizing singular configurations with zero eccentricity. The formula 
simplifies in the planar case: 


114*^11 


IXi-i -Xi| i = ,n-l 

Xn-i i = n 


where \w\ := VnP, for a given w € <C. 


g) Reflections are not well described in the framework of the Dep-reduction: Com¬ 
pare, e.g., [29, Section 4.4]. Instead, in the framework of the P-reduction, the 
transformation 

(ej) i-e,2kTr-d) 

corresponds to changing the sign of the second component of any and any 
Therefore, any of the points 

(0,^) = (0,fc7r) fcGZ"-i 

is an equilibrium point for the Hamiltonian, corresponding to a co-planar configu¬ 
ration. Compare Proposition 2.17. 
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